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^ : Abstract 

We consider type OA matrix model in the presence of spacelike D brane which is 

localized in matter direction at any arbitrary point. In string theory, the boundary 

state which in matrix model corresponds to the Laplace transform of the macro- 

_ , scopic loop operator, is known to obey the operator constraints corresponding to 

D , open string boundary condition. When we analyze MQM as well as the respective 

(~| , 

collective field theory and compare it with dual string theory it appears that consis- 
tency of the theory requires a condition equivalent to a constraint on the matter part 
^ ! that needed to be imposed in the matrix model. We identified this condition and 

Q^ , observed that this has only effect into constraining the macroscopic loop operator 

so that it projects the Hilbert space generated by the operator to its physical sector 
at the point of insertion while keeping the bulk matrix model remains unaffected, 
thereby describing a situation parallel to string theory. We analyzed the theory with 
uncompactified time and have shown explicitly that the matrix model predictions 
are in good agreement with the relevant string theory. Next we considered the the- 
^ . ory with compactified time, analyzed MQM on a circle in the presence of D brane. 

■ We evaluated the partition function along with the constrained macroscopic loop 

operator in the grand canonical ensemble and showed the free energy corresponds 
to that of a deformed Fermi surface. We have compared the matrix model features 
with that of the relevant string theory. We have also shown that the path integral 
in the presence of D brane can be expressed as the Fredholm determinant. We have 
studied the fermionic scattering in a semiclassical regime. Finally we considered 
the compactified theory in the presence of the D brane with tachyonic background. 
From the collective field theory analysis we have predicted the right structure of 
the theory in the presence of D brane. We evaluated the free energy in the grand 
canonical ensemble. We have shown the integrable structure of the respective par- 
tition function and it corresponds to the tau function of Toda hierarchy. We have 
also analyzed the dispersionless limit. 
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1 Introduction 



The two dimensional string theory (see e.g. [1], [2], [3] for reviews) is a very instructive 
model when we would like to understand the nature of string theory as a complete theory 
of quantum gravity. This theory has a powerful dual description of c = 1 matrix model 
defined by the simple quantum mechanics of a Hermitian matrix $ with the inverse 
harmonic oscillator potential U{^) — — after the double scaling limit. Matrix model is 
successfully used to describe 2D string theory in the simplest hnear dilaton background 
as well as to incorporate perturbations. 

In last decade the c — 1 matrix quantum mechanics has received lots of attention 
because of its new interpretation as the decoupled world volume theory of unstable DO- 
branes [4-6]. The matrix model dual to type string theories were also proposed in [8,9]. 
In particular, the matrix model dual of the two dimensional type string gives a non- 
perturbatively well-defined formulation. For example, the type OB model is defined by 
the hermitian matrix model with two Fermi surfaces. The type OB matrix quantum me- 
chanics (MQM) describes open string tachyons living on the unstable DO-brancs, whereas 
the type OA MQM describes tachyonic open strings stretched between stable DO- and 
anti-DO-brancs. Upon compactification on Euclidean time, these two matrix models are 
conjectured to be T-dual to each other. The exact agreement in free energy was found 
in [9]. Matrix model dual of type string in the flux background was explored in [29,30]. 
However, unlike c = 1 matrix model which can be derived from discretizing the Polyakov 
action on the string world sheet, such a derivation is not known for type matrix models. 

An attempt was made in [13] to obtain the exact form of the macroscopic loop operator 
in Type string theory. If we consider the bosonic string partition function 

/ D(I)DX exp - [ (fz[ ^{dXdX + 9090) + Qi?0 + ne^^^] - I [^ + Hee^'^ j 
J \_ J An JdT, 2n 

(1.1) 

the macroscopic loop operator inserts an operator 



W{t, I) ~ 6 i^J^^ e^-lj- 5{X'' - t), (1.2) 

within the path integral [10, 11]. 

{W{1)) = Z(l) = J DXD(l)D[ghost] s(^J^^e''' -iy6{X°-t) f{x, (J)) Z{(l){a),X (a), [ghost]), 

(1.3) 

where f is some wave function for matter ghost and Liouville. The physical meaning of 
this operator in two dimensional string theory is the presence of a 'Euclidean D-brane' 
localized in the time direction. To be more precise after we take the Laplace transforma- 
tion / d0e~'*^^'^, we get a D-brane with the Neumann boundary condition in the Liouville 
direction and the Dirichlet one in the time direction 

r dl 

J -e-^'Bi w,os{t,l) ~ \B^fzzt){i^b))^®\D)xo. (1-4) 
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when we impose the condition that boundary Liouville term is zero. 



dn<f> + f^Be'"^ = 



(1.5) 



Where (9„0 denotes the Liouville momentum normal to boundary while along the bound- 
ary we have dfX" = 0. Now consider 2D superstring action obtained from extending the 
bosonic fields to their superspace and expanding the 2D superspace action in terms of the 
component field, 



We can also consider the macroscopic loop operator which is the superspace analogue of 
Wbos{t, I), inserts the boundary condition on the fermionic coordinate iplz) = rjiplz) where 
r] — ±1 describes the RR and NS NS sector. Laplace transform of this operator inside 
the string path integral describes the boundary states NS NS and RR sector. However 
depending on helicity, in each sector we have two types of boundary states e = ± so 
that we have four types of macroscopic loop operator given by W^g, W^g, W^. The 
parameter /ib corresponds to the boundary cosmological constant in the boundary state. 
Indeed we can show this relation [33, 34] by computing one point function on the brane or 
equally annulus amplitude as shown in [12]. For c—1 matrix model the expression of these 
operators were obtained and its equivalence to string theory is verified in [10-12]. Author 
of [13] obtained the expressions of macroscopic loop operator in Type OB matrix model and 
also for NS sector of Type OA matrix model which was verified by calculating the one point 
function. Now once we understand the duality between noncritical string theory in the 
linear dilaton background and Matrix model, its natural to ask whether we can understand 
the string theory with nontrivial background which has an obvious realization in matrix 
model by adding perturbations which survive in the double scaling limit. There are two 
ways to change the background of string theory: either to consider strings propagating in 
a non-trivial target space or to introduce the perturbations . In the first case one arrives at 
a complicated sigma-model. Not many examples are known when such a model turns out 
to be solvable. Besides, it is extremely difficult to construct a matrix model realization of 
a general sigma-model since not much known about matrix operators explicitly perturbing 
the metric of the target space. Thus, we lose the possibility to use the powerful matrix 
model machinery to tackle our problems. On the other hand, following the second way, 
we find that the integrability of the theory in the trivial background is preserved by the 
perturbations. Also when we study the theory in a nontrivial background in most of the 
cases the target space metric of such backgrounds is curved and often it incorporates 
the black hole singularities. In the superstring theories, the supersymmetry allows for 
some interesting nontrivial solutions which are stable and exact. But the string theory 
on such backgrounds is usually an extremely complicated sigma-model, very difficult even 
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to formulate it explicitly, not to mention studying quantitatively its dynamics. The two- 
dimensional bosonic string theory as well as Type theory are the rare cases of sigma- 
model where such a dynamics is integrable, at least for some particular backgrounds, 
including the dilatonic black hole background. A physically transparent way to study the 
perturbative (one loop) string theory around such a background is provided by the CFT 
approach. However once we try to understand higher loops or multipoint correlators, we 
have to address ourselves to the matrix model approach to the 2D string theory . 

The 2D string theory has been constructed as the collective field theory [25], [28], in 
which the only excitation, the massless tachyon, was related to the eigenvalue density of 
the matrix field. Now consequence of the deformation in eigenvalue density correspond- 
ing to deformation in string background at classical limit was studied in [26]. C — 1 
string theory perturbed by tachyonic mode studied in [32]. Vortex perturbation and its 
equivalence to sine -Liouville theory was studied in [36], [37]. Its shown that partition 
function is integrable and have Toda structure. Toda structure and Lax formalism in the 
context of matrix model described in [31]. Many more works in this direction was done 
in [15-19,21-23]. 

Now its an interesting question to ask that can we study this nontrivial background in 
the context of dual matrix model in the presence of D brane which are just the Laplace 
transform of the macroscopic loop operator as wc discussed. Open close duality predicts 
that partition function must have integrable structure. So we consider Type OA MQM 
with simplest macroscopic loop operator, which is the operator in NS sector(as prescribed 
in [13]), which is localized in time direction and with it we show that partition function 
indeed have an integrable structure. We obtain the string equation. 

The plan of this work is in section 2 we are going to consider basic Type OA MQM 
in the presence of D brane, with uncompactificd time. We are going to introduce a no 
leakage condition to matrix model which is equivalent to some constraint to boundary 
state of string theory. We are going to explain its origin as well as its string theoretical 
interpretation. In section 3 we consider MQM compactified on a circle in the presence of D 
brane. From path integral approach we are going to show that the partition function can 
be expressed as Fredholm determinant. We have explicitly evaluated the thermal partition 
function and have shown that without application of this constraint the partition function 
diverge. In section 4 we have discussed scattering in semiclassical regime. In section 5 
we have considered string theory in the presence of momentum modes and have shown 
that the partition function in this background have an integrable structure if we apply 
this constraint. 
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2 Type OA MQM in the presence of the D-brane 



2.1 Type OA MQM 

Let us start with the MQM of type OA theory in two dimensions, which is the decoupled 
world volume theory of (stable) DO -brane and anti DO-branes. A spacelike DO -DO 
pair, i.e with Neumann boundary condition in Liouville direction and Dirichlet in matter 
direction, gives a macroscopic loop observable of the matrix model after Laplace trans- 
formation [9]. We are going to consider Type OA MQM in the presence of this operator 
and study the relevant physics. Here is a brief review of the Type OA MQM. In the 
background with no net DO-brane charges, the matrix model has U{N) x U{N) gauge 
symmetry. This is the case we are going to consider. We have the U (N) x U{N) gauge 
field Aq and bifundamental tachyon $, 

The action is 

/ dtTr IiDoMYDoM + -^M^m] , (2.3) 
J I 2a' J 

where 

DoM = doM + iAM - iMA. (2.4) 
As M is a complex matrix so we denote M by Z, — Z 

DoZ = doZ + iAZ -iZA ; {DoM^ = D„Z = d^Z + iAZ - iZA. (2.5) 

Now define 

1 _ 1 

Z±^ ^=Z ±DoZ ; Z±^^=Z±DoZ. (2.6) 

The Type OA matrix model action in terms of the light cone variable 

= y dtTr 'Z+DaZ- + Z+{DaZ) + \(Z-Z+ + Z+Z_) . (2.7) 

The gauge field A acts as a lagrange multiplier which projects the theory onto singlet wave 
functions. Its shown in [9, 15] that type OA MQM when projected to singlet sector can 
be represented by non-relativistic free fermions moving in a two dimensional upside-down 
harmonic oscillator potential 

H = \{pl+fy)-^{i' + f). (2.8) 

The theory has different independent sectors labeled by net DO-brane charge q, which is 
the same as the angular momentum J = xpy — ypx [9]. Here we will consider the case 
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where there is no net DO -brane charge, namely the J = sector. Now with z, z = x ±iy; 
and hght cone variable are as defined in (2.6 ) we have the hamiltonian 



2i 



Ho = --(z+z^ + z+z^ - 



. 2a' 

l^^^[z±^ + z±^ + l]. 



/2a'' '^dz± dz± 
The commutation relation satisfied by these operators 



[z+,z- 



[z+, z. 



'2a! 



so that 



d 
dz^ 



_d_ 

dz4 



We have Schrodinger equation 
.d 



_d_ 
dz4 



_d_ 

dzA 



(2.9) 



+ l]^{z+,z^,t). 



(2.10) 
(2.11) 

(2.12) 



Note, here we have absorbed the Vandermonde determinant in the wave function so that 
the wave function in (2.12) describes a fermion. 



2.2 Type OA MQM in the presence of the D brane : The con- 
straint from dual string theory 

Consider the type OA matrix model in the presence of D brane which arises when we insert 
an operator dtw{t)5{t-to) ^-^q matrix model path integral where W{t) is the Laplace 
transform of the macroscopic loop operator( [13], [18]). In the dual two dimensional type 
OA theory this means that there is one Euclidean DO-DO brane is locahzed at time t^. 
The branes extends along the Liouville direction after the Laplace transformation. The 
macroscopic operators can be divided into NSNS and RR sector part such that they cor- 
respond to the NSNS and RR sector part of the D-brane boundary state. Moreover, since 
we know that there are two types of (FZZT-like) boundary states |-B(e)) according to 
the spin structures there should be two macroscopic operators W^*^^-* with e = ± in each 
sector. First consider the simplest expression of macroscopic loop operator which is the 
one in NS NS sector as prescribed in [13] and expressed as e~^*^^*^(*°). Now, consider the 
Laplace transform of the operator 



^-^p^-iMm ^ _Tr log{l + (2.13) 



ZZ 

-Tr log(l + —) (2.14) 
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zz . 



= -^l„g(l+(£±±ir)(^±Iz)) 

E . Zj^Zj. + Z—Z— + Zj.Z— + Z—Zj.s^ 
l0g(l + -2 ) 

= W{z+,z+,z_,z_). (2.15) 

(Here implies sum over the eigenvalues ). Now the macroscopic loop operator for NS~ 
sector can be expressed as 

w^-y: iog(i + '^^^ ^ \ ^ ^-^^ ) 

= _ ^|gizl)![ ^+^+ +^-^- +^+^- +^-^+ ]"}. (2.16) 

n=l /^B 

The path integral over the Euchdean time in the presence of D brane is expressed as 

j XldZ+dZ_dZ+dZ_dAdAe- ^ (2.17) 

apparently impUes a shift^ ^ in the free single fermion hamiltonian (iH — )> (iH + W5{t — to) 
(in Euclidean time). Note here we arc going to consider complete quantum theory along 
with the macroscopic loop operator W. However before proceeding note that the operator 
dtw5(t-to) -^j-^icj-^ jg localized at t = to, in general breaks the time translation symmetry 
of MQM action. So as far as matrix model action in the presence of brane is concerned 
(as considered in [13]) this is describing leakage of energy exactly at tg. Now to be more 
precise consider MQM path integral in the presence of an operator e^^^°\ for which under 
any infinitesimal variation in time t — )■ t + e{t), Ward identity implies 



5(t - io)5(e^(*°)) = -dt{H,(t) e^(*°)) 

=> ( 5W(to) e^(*°) ) + lim /*°^' dt{Ho(t) e^(*°) ), (2.18) 

e^O Jto-e 

where 5W is the variation of W due to infinitesimal time translation at fixed time to- 
The above identity arises when we consider the first order(in e) variation. However e^^^°^ 
is a coherent source of the operator W{to) and on expansion generates infinitely many 

^Note that as W{z+, z+,z-,z-) in any sector expressed in light cone variable, so does not involve any 
derivative. Hence wc can just add it to hamiltonian or lagrangian as a potential localized in to 

•^Note, when we are adding the term WS{t~to) in the expression of the hamiltonian from the operator 
QW{to) jg supposed to add in the hamiltonian the terms like [/3 J^"^^ dtH, J^"^^ dt'Wd{t' —to)] + ....higher 
commutators = e~^^*''^[/3 J^"^^ dtHje^^*"^ + ...0{pef'.. . However as around the 5 — function, e can be 
made arbitrarily small i.e e << so upto quasiclassical limit one can just put these terms to zero while 
in the classical limit these terms are trivially zero. 

"^Here we have used the fact that equation of motion in the presence of an insertion 6^^*°^ is satisfied 
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source W in the path integral. So in principle we sum up the contribution from every 
order of e this gives rise an operator in the path integral 

^g[eatifc,(t)+e5W(t)(5(t-to)+{Ty(t<,)+e[ati?o(t),W^(to)]+---higher commutators}]^ 

(with proper time ordering of operators). However integrating the argument of expo- 
nential over an infinitesimally small interval to — e to to + e,time translation invariance 
implies 

/*°^' dt dt{ Hoit) ) + 5{ Wito) ) = 0. (2.19) 

Jto — t 

(one can verify the commutator terms in this integration will give zero because of 
time ordering) which is an operator constraint. Note the term VF(io) which is completely 
localized at a point to essentially creates the effect of boundary in the matrix model 
action which is defined on infinite real line in the time direction. So any variation of the 
expectation value of the hamiltonian (from (2.18)) exactly at to due to the interaction 
[Ho, W5{t — to)] with the source W{to) is the signal of leakage of MQM hamiltonian Ho 
exactly at to- So we conclude that the time translation invariance of the path integral 
implies^ that in an infinitesimal small interval around to we have ( Ho{t) )|t°lc+^(W^(^o) ) ~ 
0, where 5W{to) = [Jl°-e[Ho, W{t)6{t-to)] is the variation of W(to) due to time translation 
to to + e(to), which indicates the leakage of energy of the fermionic system. However 
the inclusion of an operator W{to) in MQM action has an interpretation in dual string 
theory is to create a boundary to string world sheet by insertion of a macroscopic loop 
localized at X° = ito or presence of a boundary state in closed string channel. So the 
above phenomenon in string theory implies that closed string hamiltonian is undergoing 
a leakage while being scattered from the boundary state localized at X = X° = ito ! This 
means that energy from the bulk is fiowing out across the boundary or in other words 
the bulk hamiltonian is not conserved in the presence of boundary! This effect can be 
visualized in matrix model from the consideration of collective field theory. 

Note the path integral with the operator e^*^*°-* does have an interpretation that free 
fermionic state are getting scattered from an operator e^^*°\ The free fermionic states 
after being scattered becomes permanently changed due to the action of an operator 
which is the function of the leakage factor / dt[SW{t)]6{t — to) at t > to- Clearly the 
scattered state will differ from the incoming state with a term which is function of 5W{to). 
As the term SW{to) is not present in the effective hamiltonian ^ or cannot arise in by 
the time evolution of (Ho) w.r.t the complete hamiltonian so it describes a leakage. 
These deformed states although evolve according to the free fermionic hamiltonian Ho 
but they can be considered as the superposition of the states which are stationary w.r.t 

^Note the effect of leakage is observed within an interval tg — e to to + e. Once we move slightly away 
from to system will evolve according to the conserved hamiltonian Ho 

6which is given by WS{t - to) + e-'^f*") [ /3 j'"^' dtff ]e^(*°) ~ pHo + W5{t - to) as we mentioned 
"^as W{to) — >■ W{to) + SW{to) is an instantaneous process 
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a hamiltonian deformed from Ho where the deformation is caused by the leakage as 
we discussed. Now in collective field theory the fluctuation of collective field from its 
static value gives a field which corresponds to 2D spacetime tachyon. The action for 
this fiuctuation gives the propagator of a 2D massless scalar. Now in the presence of 
the operator e^^^°\ the wave function of this scalar field above t > to although evolve 
according to the hamiltonian of a 2D massless scalar but will be deformed from the one at 
t < to hy the action of an operator which is function of the leakage factor 6W{to). As this 
operator is not present in the path integral so the deformation of the wave function above 
t > to must show up as the modification of the propagator from that of a 2D massless 
scalar ! This can be easily seen from the canonical quantization of 2D massless field and 
considering the deformation of the wave function. 

So we see that although the time translation invariance is maintained by extending 
Ward's principle to every order but its not giving the right string theory picture! This is 
because the closed strings which are getting scattered from D brane the scattered states 
remain the same on shell states w.r.t the hamiltonian same as that for incoming one! So 
it appears that we must need to impose a constraint in matrix model side in order to 
extract right string theory from it. Lets briefly go through the string theory scenario and 
try to understand string theoretical origin of such constraint. 

Note as far as the open string world sheet is concerned Dirichlet boundary condition 
which fixes the matter coordinate X = X° = ito at the boundary implies nonconservation 
of the momentum associated with the matter direction. However the open string action 
S'open and the open string path integral Zopen remains invariant under an infinitesimal 
variation of X which is ensured from the boundary condition 



where 5x implies infinitesimal variation in X at every point of the world sheet. Also the 
conservation of the string hamiltonian is associated with the boundary condition 



which ensures there is no leakage of energy at the boundary ^ [39, 40] . Moreover we 
have the constraint from the current algebra and superpartners of all the above condi- 
tions which ensures the conservation of the symmetry generators. Now the string path 
integral with a Laplace transformed macroscopic loop (which creates the boundary local- 
ized at X — X° with the imposed wave function giving right string one point function ) 

^To explain a bit more, in the presence of the D brane we know the energy of the incoming state is 
(associated with closed string )is not same as that of the outgoing state in the direction with Dirichlet 
boundary condition as D brane act as a source. However (2.20) implies we can consider the incoming 
and the outgoing state as the separate conserved system evolve according to same conserved hamiltonian 
(but different state) with none suffering any leakage at the boundary 



5X{X") = ^ ^x^-open = 



<^x[^open] — 0, 



(2.20) 




(2.21) 
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must obey the conditions (2.20,2.21) where the condition (2.20) follows from the function 
^(-^boundary — -^") present in the macroscopic loop functional. Hence these conditions must 
show up in the dual matrix model with a D brane. To state more precisely the 2D path 
integral on a manifold with the macroscopic loop (for the bosonic case which is given by 
(1.3) ) corresponds to a physical state where the respective wave function satisfy WdW 
equation [10]. WdW equation implies invariance of the wave function under the action 
of the generator of T(worldsheet time) translation. However the conservation of these 
symmetry generators follows from these conditions and hence the respective state must 
have information about it. So to gain the insight about what these constraints correspond 
in the matrix model let us look at the closed string channel and express the constraints 
in terms of the boundary state. Now in the minisuperspace approximation only the zero 
mode part of the constraints will be relevant and can be expressed as 

(L„-L„)|S)=0 ; (5Xboundary|S)=0. (2.22) 

Note that both the above constraints are followed by their superpartners however as far 
as zero mode is concerned we have already applied such constraints when we classified 
macroscopic boundary state according to their spin structure [9]. The second condition 
essentially describes the zero mode condition 

{X-X'>)\B)^0 =^ \B) = 6{X - X'>)\0). (2.23) 

In closed siring channel path integral with a Laplace transformed macroscopic loop corre- 
sponds to boundary states [39] . This can be seen by expressing the path integral functional 
^(X, 0) as a sum over operators(in Minkowskian signature) Oi by state operator mapping 
where we know that these operators corresponds to Ishibashi states and the wave func- 
tion {^\Oi) gives the one point function. So 1^^) must be annihilated by the constraints 
from (2.23) [42]. So the same constraints must be imposed on the state associated with 
matrix model path integral in the presence of D brane. This is because macroscopic loop 
operator in matrix model can equivalently be expressed in terms of operators along with 
the respective wave function where each component is in one to one correspondence with 
the one in string theory side. So we conclude that the kind of leakage we discussed at the 
beginning of the section is caused due to absence of any condition equivalent to (2.22) 
and must be cured once we impose an equivalent condition to matrix model. Lets try 
to find out the constraint in matrix model. From the first condition in (2.22) along with 
the one from the zero mode part of the current algebra in Dirichlet boundary state for 
matter just emphasizes the fact that 2D boundary state or the macroscopic loop operator 
will correspond to the superposition of primary states/operators expressed in terms of 
the momentum modes(i.e no winding modes) [39] with a refiection symmetry P — )■ —P in 
matter as well as Liouville part [33,39] which is already known in the matrix model [2,28]. 
The reason we obtain Liouville one point function in exact form from matrix model, is 
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that in minisupcrspacc approximation this condition is trivially satisfied and consequently 
not going to impose any constraint in matrix model side. However the second condition in 
(2.22), (2.23) or (2.20) is not yet properly understood in the matrix model. More precisely 
in the presence of macroscopic operator the state from the path integral is represented by 
an wave function expressed as a functional of bulk d.o.f ^ = ^({X}, {(j)}...) and under 
any infinitesimal transformation X ^ X + 5X we must have 

5^ ^ [ 5XJ{X)^^0 (2.24) 

J boundary 

which ensures the bulk conformal invariance and imphcation of (2.23). Note ^ is the wave 
function ^ = ^(-'^^boundary, ^boundary) which is an eigenfunction of complete string hamil- 
tonian, representing BRST invariance. Similarly its discussed in [4, 6] that matrix model 
path integral in the presence of an operator e^^*°^ ( which arises by including a probe 
eigenvalue) is an wave function ip — ip{z±z±{to) ) which satisfied the Schrodinger equation. 
However considering the fact that this operator also creates an effect of boundary and ip 
is a functional of MQM variables ip — ip{{z±z±{t)}) we must have an condition equivalent 
to (2.24) in matrix model which ensures conservation of MQM hamiltonian in presence 
of such operator. Naturally no such constraint arises from Liouville d.o.f for the reason 
as we discussed. 

Here first we will find such constraint in the matrix model from somewhat intuitive 
way, solve it in the context of the matrix model path integral and show its consequence. 
Finally with the help of it we will come to more rigorous analogy between the string 
theory and the matrix model scenario in the next subsection. 

First note that in the matrix model the matter coordinate X is getting mapped to 
time(Minkowskian) coordinate, X — > it. So we can guess that string theory boundary 
condition must be refiected in MQM as an overall invariance of the path integral under 
time translation i.e SX = 6t with no leakage. Note when there is no leakage, under 
infinitesimal time translation t — )■ t + e(t) the variation of path integral is given by 
((5e-/^'^*'^^*-*'^^*~*°-'), where d defines the variation of the operator due to infinitesimal time 
translation at fixed time t — tg- So string theory boundary condition (2.20) must be 
refiected in the following constraint in matrix model 

^jdtwmt-to)^^ = dtwmt-to)^ ^ 0. (2.25) 

Indeed in the string theory path integral if we expand the macroscopic loop in terms of 
operators which corresponds to Ishibashi state one can verify that this is the consequence 
of Ward identity under an infinitesimal transformation X — )> X + 6X which arises on 
application of the second condition in (2.20) and we have already mentioned it in an 
alternative way in (2.24). 

In next subsection we will show that the consequence of this condition are in exact 
agreement with that of string theory. To understand the impact of this condition in MQM 
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first we need to write down the Schrodinger equation and study the Hilbert space. We 
have the time dependent Schrodinger equation for a single fermion ^ in Minkowskian time 

d 



For t away from to we have time independent Schrodinger equation 

ig-^'4^{z±z±, t) = + + l]^^f{z±z±, t) = Eij{z±z±, t). (2.27) 
The free fermion solution with energy E is 

= e-^^^e^^'^iz^z^r'-K (2.28) 



where we have chosen a' = 2 and (f)o{E) is determined from biorthogonal property (dis- 
cussed in the Appendix) and given by 

e^ME) ^ p-^ + i) (2.29) 

Now consider the commutation relation 

[Ho,z+\ = -i ; [Ho,i+] = -i, 

[Ho,z.]= I ; [i/„,t_]=i. (2.30) 

This imphes and when acts on a state \E) expressed \E) ^ \E -i) 

and Z-Z-\E) — \E -\- i). These states can be represented as \E :kni). These describe 
a different Hilbert space [20] which can be understood as their inner product with the 
state \E) either diverge or zero. These states actually can be identified with the discrete 
tachyonic states over the matrix model ground state [20,28]. 

Now to the meaning of the constraint. First we will find the expression for the con- 
straint and then solve it in the context of the matrix model path integral for the macro- 
scopic loop operator we considered or in the more complicated case to reach to the right 
expression free energy. Consider the v.e.v of the operator in single fermionic state which 
is given by 

(e^(*")) = = ((1 + !±f±±Zr!^^)^. (2.31) 



^When we consider the insertion of j3 factor, for the macroscopic loop operator we have the expression 

■gFT, however when we rewrite the Schrodinger equation in terms of the eigenvalues x,y which corresponds 
to the real and imaginary part of eigenvalue z we have the Schrodinger equation 

In the double scaling limit we take x,y ^ \/^x, y/^y [1], which gives the Schrodinger equation (2.12 ) 
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(Here we have shown the expectation vahic w.r.t the single fermionic state. This we could 
do because the theory is projected singlet sector and N fermionic state is just the direct 
product of each). So (2.25) along with (2.30) gives the constraint 

{z+z+{to) - z^z^to)) = 0. (2.32) 

As the constraint is exclusively on the variable associated with W so it must has effect in 
constraing the ffilbert space created by W at to- Note for the macroscopic loop operator 
in any other sector, in general the variation of {e^^^°^) can be expressed as 

{e^(io)sw{Q) = (e^(*°)[^+z+(i„) - z.z.iQ]) = 0. (2.33) 

As the path integral with the uncompactified time essentially describes the transition 
amplitude from the initial state \'z±z±{ti)) to the final state \z±z±{tf)) , so the condition 
(2.25) remain true for any arbitrary variation in time implies that we have 

(^;|5e'^(*°)|V^,)=.0, (2.34) 

between any two physical states which according to the Hilbert space described 

by (2.30) will be of either form \z+z+,t) or \z-Z-,t). Now recall the expression of the 
inserted macroscopic loop operator W{to) in (2.15) which is expressed in the form (2.31). 
While the Hg in the expression of W keeps the underlying state invariant the ^+-2+ and 
act on the vacuum to create the states \E = —i) , = i) or in other words the 
macroscopic loop operator (2.15) due to its 'zz field acts on the vacuum to create a left 
moving as well as right moving state. So the constraint essentially relates the two exactly 
at to by putting the following constraint inside the path integral, 

z+{to)-z.{to) = ■ z+{to) - z.{to) = 0, (2.35) 

The physical meaning of the above condition is that the quantum fluctuations of the 
variables z+{z+) and z^{z-) around its classical value appears to be identical at tg- This 
in turn implies that the wave functions of the states which are being created by the action 
of Z+Z+ and z^Z- over the ground state in z^ and z^.z^ representations respectively , 
appeared to be indistinguishable exactly at to- Note that in the above constraint (2.35 ), 
the l.h.s will not remain invariant once we move away from to- In otherwords if we express 
them in terms of the respective operators(which acts on vacuum), we see that l.h.s does 
not commute with the free hamiltonian Ho (2.9). So once we move away from the point 
t = to Viie have our physical Hilbert space described by the free hamiltonian Hg and the 
relation (2.30) so that again the left and right moving states created by W will appear 
to be distinguishable. In next few steps we will see that the condition (2.35)whcn acts 
inside the path integral it just has a meaning to project the operator W and the Hilbert 
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space generated by it to its physical sector while keeping the bulk physics unaffected! So 
with the constraint we can express the path integral somewhat schematically as 



I 



t<to 

dz+{to)dz+{to)dz_{to)dz_{to) d{z+{Q - z^toM^+iQ - ^_(iJ)e^(^+^+'^--'^°^*°) 



/ n dz+dz+dz^dz-e--^-oo'''' = / dz+dz+dz_d^_e--/'--'^*[''^-^(^±^±(*°)'^'')l. (2.36) 

t>to ■' 

Its important to note that the effect of the constraint is only to project W at its physical 
sector without imposing any boundary condition to original lagrangian which happens 
due io 5{t — to) factor as we explained. Projection implies we can get same path integral 
expression or the transition amplitude by expressing W either in z+z+ or in z-Z- mode 
which happens due to the fact that wave function associated with the either mode appears 
to be same at to- Its also important to note that when we consider the original expression 
of W (2.15), the quantum fluctuations of the varibles z+{z-) and z+(z_) are constrained 
by (2.35) and in that sense when we consider the complete Hilbert space as described 
by (2.30). -2+(-2-) and z+{z^) are not ordinary operators. However the implication of 
(2.36) is that we can replace the theory with the expression of W (as in (2.15) ) by the 
projected one W — )■ VFproj = W{z±z±, Ho) and in the theory with Wproj these operators 
act as ordinary operators and one can evaluate the partition function in a formal way in 
MQM by using W^proj instead of W to give the right transition amplitude as in the original 
theory with the constraint(2.35). 

Now we will show that how the constraint is solved in the context of the path in- 
tegral leads to the right expression of the partition function. Note in (2.34) \ipi,ti) can 
be expressed in the free fcrmionic Hilbert space \E) while the constraint (2.25) implies 
the macroscopic loop operator which has expression W{z+z+(to), Ho) or W{z_z_(to), Ho) 
leads to the most general expression of the final state \ipf) is Y^Co-e\E) + J2Cn-E\E + ni) 
or J2 Co;e\E) + J2 C-n;E\E — Tii) . As the tachyonic states have imaginary energy so the out 
state with indefinite number of tachyons will contribute to the path integral. So (2.34) 



above identity can arise by inserting 6l(zj^{to) — z-{to))5{z+{to)—z-{to)) in the path integral, operating 
Z-Z- = al^gf^ on the inner product {z+z+{to + dto)\z+z+{to))- Its important to note that due to 
S{t — to), inside the path integral we can make the time interval around to arbitrarily small i.e dt « ^ 

so that around to, e ■'*<'-« will just be identity and hence the effect the constraint is only to modify 
W without affecting the MQM lagrangian L and hence no boundary condition will be imposed on the 
variables in original lagrangian a,t t = to- Finally integrating over z_(to),'z-(io) where the S — function 
gives J-jZ- = 2+ , 2:+ , leads to the above expression of W. As W becomes independent of Z- , Z- so we 
can just express the integral as a continuous integral. Also note in the expression of Ho),the 
action of Ho on any intermediate state in the path integral gives the constant and it can be alternatively 
given by ^. So it will not be affected by the constraint. 

^^note that here path integral is expressed in a schematic way omitting the angular factors. One can 
verify that even the inclusion of angular factor will not change the picture 
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now can be expressed as a trivial identity 

oo oo 

^ c^n-AE-ni] /|exp[iy(z+z+, {to))\\E) = ^ Cn-AE+ni- /|exp[iy(z_z_, {to))]\E). 

n=0 n=0 

(2.37) 

As ip{z+z+, E) and are the same wave function in different representation 

so we see in the above l.h.s and r.h.s are the same expression, expressed in different 
representation. Now before coming to the string theoretical interpretation of all the 
matrix model events first we will derive the expression of the wave function (physical) 
for t > to- Integrating (2.26) over the infinitesimally small interval around t — to and 
following the constraint from (2.35,2.36) we have 

i-^{z±,Z±,t)\to+e - i-^{z±,Z±,t)\to-e = iWproj*(^±, to) (2.38) 

I l0g(l + \ -)'^{Z±,Z±, to) 

(2.39) 

So exactly ai t — to the wave function is given by 

limj/ ^iz±,z±,to) = (1 - Wi^±z±,Ho))^o{z±,z±,to) (2.40) 

So for the macroscopic loop operator in the NS NS sector we have the physical wave 
function for t >to 

^A-^^z^,to) = [1 - log(l - '^'^'^^^"^~'^^" )]^o{ziZ^,t) (2.41) 

Now we need to find the wave function at t > tg- In order to do so first note that for 
t > to wave function evolves according to the free hamiltonian Ho- So at the first sight 
it appears that the wave function at t > to is given by the one obtained from the time 
evolution from '^>{z±z±,to). It is given by 

^>{z^z^,t) = [1 - log(l - ?^^^^^^^)]^o{ziZ^,t) (2.42) 

Now although exactly at t = to, W is expressed in the form W{z±z±, Ho) but once we 
move above to, W can in principle be expressed in terms of both ^+-2+ and 'z_z_ as both 
of them are related by the constraint (2.35) at to- More precisely these states are given 
by 

IV'/) = Ec™n(i?)(t+5+)"(t-^-)"|i?), (2.43) 
So following the previous discussion \ipf) is expected to be given by 

{l-W{^+z+{t),%+z4t)Ho)),^o{z±,z±,t). (2.44) 

We can express these states as the one time evoluted from the wave function from t = tg. 
Although it appears that exactly at t = to, \ipf) and |^>) are of similar structure but they 
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are not same as time evolution property of ^+ and 5+ are different and once we move away 
from to we liave our original Hilbert space. Finally note that at t = to fermionic states 
are being converted from ipg to so the change of fermion number N must be measured 
from the variation of the transition amplitude at t — to i.e in terms of the macroscopic 
loop operator described in (2.15) we have 

d , , , M , / XV ,dW . , 

Where ti — >■ — oo and tf — >■ oo. The average {^)\t=t„ — described above is w.r.t the 
path integral. So in order to have the fermion number unchanged we must need to impose 
the condition 

(^>l«. = 0. (2.46) 

Note that this is just a condition parallel to (2.25) and have an effect to project W as 
described above. In the next part of this section we are going to see that all these matrix 
model events are exactly in one to one correspondence to the string theory. 

2.3 String theoretical interpretation 

In this section we will see that the constraint (2.25,2.35) we imposed leads to right 
string theoretical result. First The boundary state of 2D spacelike brane is expressed as 
|-B(SupcrFZZT)(/^B))0®|-D)x" where |-D)xo for NS NS sector [39] is expressed as A/" dPe^^^° 
( J\f is the normalization factor) which in terms of the vertex operator can be written as 



\B) 



/oo 
(iP^[e'^-(^-^°)] |0) + descendants. (2.47) 
-oo 

The string cndpoints are localized at X = X° and the matter part of the boundary 
state form the representation of 6{X — X°)\0), satisfying (2.22). Boundary state of Super 
Liouville NS NS sector is given by 

(2.48) 

where \vp^) primary macroscopic state associated with a vertex e*^^+*^'^''^, M is the mass 
of intermediate propagating mode. 

The boundary state is the direct product of the matter and Liouville part along with the 
ghost factors. Setting Pm = Pi = P we have the primary part of the boundary state 
without ghost excitation mode, which is the superposition of the tachyonic field can be 
expressed in terms of the operators from the state operator mapping as 

J ciP[/(P)[e^^(^-^°) + e-^(^-^°)]e(*+^^)^ (2.49) 

Now in the 2nd quantized matrix model we can express the macroscopic loop operator as 

W{l,t) = J e^~'ij^ = J dTe^'"''^'^drV{l,t), (2.50) 
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where r is the time of flight coordinate obtained from reparameterization 1.2 = ; ~ 
cos^/iT , ~ dTri{l, t) and rj is the massless bosonic field {d^ — d'^)ri(t, r) = which cor- 
responds to the tachyon in the string theory at asymptotic r [2,25,28]. rj corresponds to 
the fluctuation of the collective field = (po + drTj. rj satisfies Dirichlet boundary condition 
in T direction [11], [28]. However note that when associated with the Laplace transformed 
macroscopic loop operator operator W which correspond to D brane boundary state, 77 is 
no longer an ordinary state but it corresponds to an Ishibashi state. Here we show that 
the constraint we imposed (2.25,2.35) leads to the right matter one point function from 
77 what is expected from string theory. Consider 



r^{r, t) = /" ^f){p){a{p)e-'P' + b{p)e'P')sm{pr). (2.51) 

J —oa p 



In order to find the exact t dependence consider the fact that the implication of (2.46) in 
the context of collective field theory is 



0, (2.52) 



where as before 5t implies variation of the collective field due to infinitesimal variation 
in t, at fixed t. 

So (2.52) implies in (2.51) we have 

^^gipto . i) = e-ip*o (2.53) 

So (2.51) is expressed as ■q{T,t) = /f^ ^fi(j)^[Q-w{t-to) _^ e*^(*~*°))sinpT. Now integrating 
over T we have W{p, t) = ^^pKip{ijLl)fj{p){e''^P^*''*°^+e^P^*~^°^) where Kip{^l) is the Bessel's 
function which is the macroscopic wave function satisfying WdW equation [13]. This on 
Laplace transform / ^e~^s' can be expressed as 

W{p,t)^U{p)f,{p,t) ; U{p) ^ ""^hiZ) ^^'^^^ 

where we have 

//g = 2sinh^(7rs)|/i| (e ■ sign(/i) < 0), /i^ = 2 cosh^(7rs)|/i| (e ■ sign(/i) > 0). (2.55) 

Note (2.25) and (2.52)are just the same constraint expressed in the Irst and 2nd quantized 
formalism. Given the form of 77, on Laplace transform of (2.50) and on inverse Fourier 
transform of (2.54) we can express W as 

W,{T,t) = t/(p)r;(p)[e^^(*-*°)+e-^^(*-*°)] = U{p)[drr],{T+{t-Q)+drVp{r-{t-Q)] (2.56) 

where the suffix p implies pth component. We know 77 corresponds to the space time tachy- 
onic field in the asymptotic r region which in the string theory side describes asymptotics 
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of Liouville field which describes vanishing Liouville walL So we see that each com- 
ponent with momentum p in the expansion of W describes the tachyonic operator fj{p) 
dressed with Liouville and matter wave function (2.47,2.48) where W is symmetric under 
P — > —P. So these are just in one to one correspondence with the state obtained in the 
expansion of the boundary state in (2.49). So (2.56) is just the same as (2.49) when we 
identify the physical d.o.f. More interestingly if we try to extrapolate rj in the region 
Pm Pl with the same r dependence as in (2.51), from (2.52,2.51,2.56) we have 

/fjl poo poo 

je-^Biw^l^ = [ dPiU(Pi)rji(Pi)] x[J ^ dP^e^^-(*-*°)] 

roo 

= [ / dPiU{Pi)7ji{Pi)] X 6{t - Q, (2.57) 
Jo 

where rji implies Liouville part of rj. This has a very similar structure with that of 2D 
boundary state 

\B) = I^Liouville) ® l^matter) = / dPiU {Pi)\P) ® 8{X - X'')|0). (2.58) 

Jo 

However this is a mere extrapolation without any physical justification. Finally we see the 
effect of the constraint (2.25) is only to give the correct form of the matter wave function 
while keeping the Liouville wave function unchanged as we discussed in section 2.2. 

Finally lets explain the meaning of the projection W — > VFproj as described in (2. 36), (2. 37). 
Note both Liouville as well as the matter part of the boundary state is symmetric under 
interchange of P — > —P (2.48). So projecting the theory to either right moving or left 
moving part i.e in the boundary state keeping either the left moving or the right mov- 
ing part while expressing the superposition of left and right moving state as only left 
or only right moving state (i.e flipping the sign of p where necessary)in the incoming or 
the outgoing sector yield the same transition amplitude as the original one. Overall time 
translation invariance of the projected theory follows from the consideration of either left 
or the right moving part as the image of the other. 

Lets come to a precise analogy between MQM, 2nd quantized matrix model or collec- 
tive field theory and string theory, which arises as a consequence of the constraint. We 
know the macroscopic loop operator can be expanded in terms of microscopic operators. 
To see that in dual matrix model consider VFproj(^±-S±) Ho),a,civag on ground state 

log(l + + - ^^<^ )\y,) ^ E«^n(/.)[^]-[^]-|/.) (2.59) 

Essentially the above which is described in Minkowskian time corresponds to a coherent 
state. Also the above is an expansion in terms of microscopic operator. This is evident 
from the correspondence between the operator (^±-2±)" and the vertex operator for discrete 
tachyonic states. [28]. From the previous discussion it follows that when wc consider the 
space time wave function for each component in W, we found an effect of boundary 
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exactly at X = X°(or, t = to). More precisely at this point the state associated with 
each component of boundary state in (2.49) or the states appear on expansion of W 
(2.56) in collective field theory, expressed with an wave function which has only Liouville 
part. So the left and right moving state appears to be identical at X — X°{ort = to). 
Now in order to understand this effect in MQM note that the macroscopic loop operator 
W{z+z+,Z-Z-, Ho) acting on ground state, essentially create a left moving as well as right 
moving state (2.59). The constraint (2.35) inside the path integral have an interpretation 
in MQM that a,t t — to we have 

[t+i+]"||/.)^[t_i_]"|/.). (2.60) 

i.e the left and right moving states which arc created by the action of W on MQM ground 
state appear to be exactly identical at to (but different when we move away from to). This 
leads to the projection W — )■ Wproj = W{z±z±{to), Ho) in the path integral (2.36). More 
explicitly the created states obey ^jJ+{z+z+{to)^, E = ni) = il)-{z-Z-{to)]E — —ni). This 
is supported from (2.35) and give the justification for (2.37). Again lets emphasize on the 
fact that these constraints no way affects the free fermionic states as well as the discrete 
tachyonic states of the original Type OA matrix model. This distinction just corresponds 
to that of the closed string states associated with the boundary states and the free closed 
string states in the bulk. 

Now once move away from to we are back to our original scenario described by (2.30) 
with the distinguishable states described by (2.41,2.44). An important point to note here 
is that when a bulk operator approaches to the boundary we encounter a singularity [40] . 
Similarly in MQM we have the ordering ambiguity of the operators arises when we move 
from the region t — to to t to on time evolution of (2.40) to (2.44). If we consider the 
entire space of N fermions N ^ oo this leads to singularity. Now in order to find the 
exact coefficient of superposition of the states from (2.41) and (2.43) above to we must 
remember the fact that these states just show up as indistinguishable at to while away 
from to they are distinct. So the coherent state above to must be given by (2.44). Prom 
the discussion of [2,11,28] these operators corresponds to special tachyonic states and 
higher Virasoro primaries. 

Now in this context we need to identify the open string operator or the boundary 
operator. Note that exactly ait = to for the states/operators from the string theory (2. 47), 
collective field theory (2.56) or from the MQM (2.60) (when we set the origin of time at 
to), they do not show any explicit time dependence exactly at to. These operators can 
either be viewed as the one obtained from the one at t ^ to by time evolution as we 
discussed throughout or the one without any explicit time dependent part. The time 
independent one must correspond to the states/ operators which are extended in Liouville 
direction but localized in matter direction. Hence they should be identified with the open 
string operators. 
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Finally let us briefly tell about matrix model string theory correspondence of the 
complete scenario we just obtained. Theory describes free fermionic state as well as 
coherent state where the coherent state is strongly locahzed at t — to- The closed string 
state turned into a coherent state at to as described by (2.39), which follows the view 
of [5,44]. The transition amphtude can be read from collective field theory. 

As the collective field correspond to a spacetime tachyon so from (2.56) we obtain the 
closed string emission/absorption amplitude, which is giving the one point function in 
mini superspace approximation 

{W V{P)) = C/(P)e-^^*°- (2.61) 

where V(P) is the operator representing tachytonic vertex and the above relation can be 
viewed on evaluation of the above in 2nd quantized field theory [2]. Now the state in (2.41) 
while show up as time dependent state w.r.t the free hamiltonian Hg, they are stable w.r,t 
an effective hamiltonian — Hot^ . The effect of the macroscopic loop operator 
is to cause transition of a fermion from Fermi level to above. At the double scaling limit 
/3 — > oo the transition amphtude for a single fermion is given from (2.28,2.40) by 

{^l^^{E\tf)\^l^{E,U)) ~ {'>l^o{E')\W\^l^o{E)), (2.62) 

with ti{tf) = — oo(oo). The presence of D brane will change the Fermi level fi — )■ /i'. So 
from the 2nd quantized theory we can see that the transition amplitude — )■ is in 
accordance with [4], [38]. 

{l^'\W{l,t)\i,) = e-^(^-'^')e^('^-^')*X,(^_^,)(^). (2.63) 

Note the transition amplitude is time independent which is evident from matter one point 
function (2.56) and it signifies the fact that we have stable D brane. 

3 Type OA MQM on a circle in the presence of D 
brane 

3.1 Evaluation of the free energy of Type OA matrix model on 
a circle 

We consider Type OA matrix model compactified on a circle of radius R. As considered in 
the previous section, there is no net background DO brane charge and hence it is described 
with U{N) X U{N) gauge symmetry. The partition function in terms of the hght cone 
variables is given by 

J dZ+dZ_dZ+dZ_dAdAe-^ir'''^'i'^+^^^-+^+'^^^^+^^^^^ , (3.1) 
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where DaZ = dtZ + i[AZ — Z A] and A is as given in (2.1). Now we fix the gauge dtA = 0, 
which sets A and A to their zero modes A'^^^ = X/2na' and /l*^*^) = X /2'Ka\ where in the 
T-dual theory X and X corresponds to collective coordinate of DO and anti DO brane [9]. 
As before the gauge fixing introduces the FP determinant [15] 



/ 



d6dcexp(™a,Ac) = nf"?'^-'^^''Vf"'f- (3-2) 



i<3 



{xi - Xj)R/2 J \ {xi- Xj)R/2 



where Xi and Xi are the eigenvalues of X and X respectively. Now the denominator gets 
canceled with the respective Vandermonde determinant so that the usual measure factor 
A(x)^A(x)^ is converted to the measure for unitary matrices 

Note these are the measure for unitary matrices U — U — e^. which are holonomy 

factors (we chose a' = 2). Therefore the natural variables to be integrated over are the 
"holonomies" U — e^*'^"^^"'-^, U — e^'^*"^*^"'^. Once we gauge fix A and A The partition 
function depends on the gauge field only through the global holonomy factor, given by 
the unitary matrix 

n = f'e^-/o'"''^(*)'^* ; n = f'e'-/o'"''^(*)'^*. (3.4) 

In the A — const gauge, in the path integral, the constant modes of A can be absorbed 
by redefining the fields Z_ , Z_ as 

Z_{t)^e-'^'Z_{t)e'^' ■ Z_{t)^e-'^'Z_{t)e'^\ (3.5) 

which replaces the periodic boundary condition 

Z±(27ri?) = Z±(0) ; Z±{2t^R) = Z±(0) 

by a 5C/(A^)-twisted one [14], [17] 

Z+{27rR) = Z±(0) ; Z+{2nR) = Z+{0) 
Z_{27rR) = nZ_{0)n-^ ; Z_(27ri?) = f)Z_(0)f)-\ (3.6) 

So in the constant A gauge integration with respect to the fields Z±{x), Z±{x) is Gaussian 
with the determinant of the quadratic form equal to one. Therefore it is reduced to the 
integral with respect to the initial values Z±,Z± — Z±,Z±{0) of the action evaluated 
along the classical trajectories, which satisfy the twisted periodic boundary condition 
(3.6). Therefore the canonical partition function of the matrix model can be reformulated 
as an ordinary matrix integral with respect to the hermitian matrices Z+ , Z_ ; Z^ , Z_ 
and the unitary matrices Q, Q, : 
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where we denote 

q = e^'^^. (3.8) 
Now note the above expression can be written as 

= J dZ'^dZ_dZ'^dZ_dfldfl'e'^^^'''^'+^-+^-^'+''^^^ (3.9) 

where we define Z^tlil~^ = Z'j^ ; Q£l~^Z^ = z'_^_ ; ^2' = r2r2~^.The last expression 
imphes replacing X hj X — X, both running over the infinite real line. The redefinition of 
the variables will keep the measure invariant. So by generalizing Harishchandra-Itzykson- 
Zuber integral we can write the above partition function as [14]. 



°? N 



[det (^e'"'^*^+^'^-*) det^-fc (^6-*"^+^^-*) det^-fc (e'^' 

det,-,(e-^^^+^-^-'=)], (3.10) 

where Q'j/^ which has only diagonal elements nonzero. Now we show that the grand 
canonical partition function can be written as a Fredholm determinant 



where 



Z{fx, t) = Det(l + e-^''^'^''K+K^), (3.11) 

[K+f]{z.z.) - l[dz+][dz+]dte'^'~'+'-+'-'^-'+^f{z+z+), 

[K^f]{z+z+) ^ J[dz_][dz_]e-''i'^'^+^-+'-'+^f{z_z_). (3.12) 



-^^(/ dUe^'^'^^-^^ = Const. ^^2^ where Xk and yi are eigenvalues of X and Y and A{x) A(y) are 
Vandermonde determinant given by A{x) = ^[{^k — xi) 

i<j 

^^To express the part involving O' we used the fact that in the integral (J dUe^'^^-^^ where D is 
a complex diagonal matrix with = and Y = VyV~^ where y is the eigenvalue of Y and V is the 
unitary matrix diagonalizing Y. Now we can write DY = V'dyV'~ for some other diagonalizing matrix 
V which exploits the fact that Di&g{{DYf DY) = Diag(yty) = Di3,g{Vy*yV^) which implies the above 
expression (d is eigenvalue of D).So following the formal derivation of the integral we can write 

(/ dUe^'^^^-^^ = Const. which is nonzero only when k = I. for any diagonal matrix 
D. Note, we are not summing over k and l.Denominators gets canceled with the Vandermonde determinants 
appearing from Z^, Z^, Z-, Z-. 
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K+K_f{z+z+) = j[dzJ\[dz-]dte'^^'+'-+* '^-^+) j[dz+'][dz^']e-''i^'+'^-+^-'+'^f{z+'z+'). 

(3.13) 

Note that t and t'^ denote the diagonal elements of Vl' and Vl'^^ corresponding to 'z±,z±. 
Now note when we evaluate the determinant in a diagonalizablc basis which is naturally 
given by f{z±z±), K^K_f{z±z±) will be independent of t i.e / dt will come out as an 
overall factor. So following the analysis of [14] the grand canonical partition function 
^ g-27rii/3Ar//2:^ can be expressed as the Predholm determinant (3.11) which is same as 

N 

that of c = 1 matrix model. Now following [24] we can express the partition function as 
Trexp[—2nR(3H]. The gauge field A project the theory to singlet sector so that in the 
absence of perturbation, the grand canonical partition function is given by the Fredholm 
determinant 

Zifi) = Det(l + e-2-«/3('^+^o)), (3.14) 

which must be same as (3.11). This can be interpreted as the grand canonical finite- 
temperature partition function for a system of non-interacting fermions in the inverse 
Gaussian potential. The Fredholm determinant can be computed once we know a complete 
set of cigenfunctions for the one-particle Hamiltonian Ho- Now in order to evaluate the 
free energy we need to find the density of states, it is conventional to introduce a cutoff 
A. 

There is no momentum flow through the wall zz = x"^ + y"^ = |Ap is implied by the 
condition ixpy + xpy)ip±{x,y)\(^2^y2^^j^^2) = {%+z+ -%-Z-)^±{z, z)\(j^=ia\^) = 0, which has 
a solution 

^^(A)=^^(A). (3.15) 
This condition is satisfied for a discrete set of energies En{n e Z) defined by 

0o(^n)-^nlogA + 27rn = O. (3.16) 

From (3.16) we can find the density of the energy levels in the confined system 

as derived in [14] Now we can calculate free energy ^(/x, R) = logZ(/x, R) as 

/oo _ - 

dE p{E) log 1 + e-2'^-«/'('^+-^) , (3.18) 
-oo '- -' 

with the density (3.17). Integrating by parts in and dropping out the A-dependent piece, 
we get 

■^(''' ^> = "S^ / *(^' + ^ _^ £ (3.19) 
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We close the contour of integration in the upper half plane and take the integral as a sum 
of residues. This gives for the free energy 

T^-i Yl (koiir/R-lA. (3.20) 

r-=n+i>0 

As the Predholm determinant is similar to that of c = 1 MQM so following the analysis 
of [14], [15] we can see Prom (3.20) it follows that [19] 

'^^^'^^■H^^)-^o{-^^). (3.21) 

Also its shown that the free energy can be expressed as 

J-pert(/i)fe=o} = -f /^^ log I - log f + E t^'-''cn{Rl (3.22) 

h=2 

3.2 Free energy of Type OA matrix model on a circle with D 
brane 

In this section we will consider the type OA matrix model path integral in the presence 
of a D brane and show the grand canonical partition function can be expressed as the 
Predholm determinant. We consider the brane in the NS NS sector and show that how 
to generalize the analysis for the brane in any other sector. Consider the path integral in 
the presence of the macroscopic loop operator localized at to, which is the generalization 
of (3.1). The classical action will remain periodic even in the presence of D brane so we 
can express the 

j dZ+dZ_dZ+dZ_dAdAe-^C''''''^[^+^^^-+^+^^'^-+^^^^^^ (3.23) 

The macroscopic loop operator depends on diagonal elements only, so the partition func- 
tion (3.7) can be expressed as 

7 ^ -1 - 

2N{t)= j n [dz+,][dz^,][dz+,][dz.,][dt,]det^ 

— oo 

det,fc(e-*^^-.^+^) explY, log{l + ^ ^-'^'^ ^ ^^'^'^ ^ ^''^^^ )], (3.24) 

(when we have off-diagonal t is zero , also wc arc not summing over j,k) 
where E 1^ coming from Trace and again above can be expressed as 

i 

oo 

— oo 

detjk (e-^^^-^-^+^) n(l + ^+rZ+r+^-rZ-r + ^+,^-, + ^-,^+, ^ _ ^g ^g) 
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Now in order to write the above expression we have used the fact that the classical action is 
periodic even in the presence of the macroscopic loop operator Z{2tiR) = Z{0). However 
at the quantum level there is a discontinuity of state \'ijj{27rR — e)) ^ |^(0) + e). This 
causes the absence of the vortex d.o.f. Now if f{z±z±) is the function which form the 
representation of and K_ (3.12, 3.13), the action of W on f is given by 

Wfiz^z^) = (1 + ^"^"^^"i-^"-'^ )f{-z^z^). (3.26) 

Now, the operator W does not introduce any interaction between the fermions, so no off 
diagonal terms from W. Now from (3.11) and (3.24) we can write the grand canonical 
partition Ejv e^^'^^^^Zjv as 

Z{i^) = det(l + e-'^^^^^'WK). (3.27) 

Now in order to evaluate (3.27) following (3.13) a representation of K is formed by the basis 
/(z+jZ+j) with i runs from 1 to N. Also from (3.13) it follows that f{z+z+) ~ {z+z+)"'. 
So when we evaluate the expectation value of WK in this basis, in the expression of W we 
see that {z±z±) = as on a closed contour the angular integral will vanish, in the inner 
product, the other term z+Z- + z+Z- expresses nothing but the hamiltonian of which f 
is an eigenf unction. So if ipn are the set of functions which diagonalizes K we can write 
(3.27) as 

^log(V'n|(l + e-''^^^^WK)\^l;r^), (3.28) 

n 

where 

WK+K_f{z+z+) = J [(iz_][dz_]e'(^+^-+^-^+) [dz+']e-^«(^+'^-+^-^+') 

(1 + H±:^^±^)/(W)- (3.29) 

As the expression depends on {z^z^ + z_z^) which is the expression for free hamiltonian 
Hg so comparison with (3.14), Fredholm determinant is expected to be given by 

Z(/i) = det(l + e ^). (3.30) 

This is, we are going to analyze in the next part of this section. 



3.3 Evaluation of the thermal partition function 

In this section we arc going to study type OA MQM in the presence of D brane with 
time t compactificd on a circle, evaluate and analyze the free energy. In the absence of 
the brane when we compactify string theory on a circle of radius R, in the dual MQM 
the Schrodinger equation have periodic solution i.e ip{t) = ip{t + 2ttR), which implies 



25 



E = ^. Now consider the theory with D brane which can be accomphshed by inchiding 
a macroscopic loop operator iocahzed at t = to = = 2nR (say) to the action. From 
previous discussion it follows that in the presence of the operator Schrodinger equation 
will have well defined solution only in the region < t < 27ri? when discontinuity occur 
at the respective point and we have ijj{2TrR — e) 7^ ■^(27ri? + e) in the limit e — )■ 0. This 
is consistent with the fact that the presence of a spacelike brane breaks the winding 
symmetry and apparently the theory correspond to that of an open string. At the end 
we will see how the closed string scenario arise in this picture. Now in a compact time 
we must have the condition t/j{t) — ■0(27ri? + t). So effectively we can view the theory as 
MQM on a hne of length 27ri? with two 5-potential along with the operator W (where 
one is the image of the other, superimposed) placed at its two ends. When we cross the 
boundary on either side, situation repeats i.e we can define the theory on any of the 
slices 27r(n — 1)R < t < 27rnR. So effectively we have the time dependent Schrodinger 
equation with double delta potential well as: 

[i§^ - {S{t) + 6{t-27rR)}W{^±z±{t),Ho)]^{z±z±,t) 



= Ti 

We have the discontinuity 



d _ d 

^ 

oz± oz± 



^iz±z±,t), (3.31) 



iP{2ttR + e) - 'ip{2nR - e) = Wijjo{t = 2ttR). (3.32) 

In order to evaluate the partition function we must need to know what is the right Hilbert 
space describe the wave function ip on the circle. This is because we know that the 
Hilbert space {ji?)} and ± ni)} cannot be mapped to each other. So, to answer this 
note that when we define the Schrodinger equation in the double delta potential well in 
an uncompactified direction we have one type of the solution inside the well, while the 
solution at the left and the rightside of the well differs from the same, decided by the 
discontinuity (3.32). Now compactification on a circle of length 27tR imply the outside 
region of the well is just squeezed to a point t — = 2t:R and the wave function at the 
left and right side of the double delta-well are given by the wave function at the right 
and left e-neighborhood of that point. Now in uncompactified time we had free fermion 
wave function (2.28) in the region —00 to to while from to to 00 the wave function is 



"^^note it never implies periodicity, its just similar to the situation of an open string in 2D with Dirichlet 
boundary condition in compact direction and identification of the matter direction with t. It winds along 
the circle m times although ends are not identified. The open string which wraps m times a circle of 
length 2-kR' with 2'KraR' = 2ttR, we can define same theory on either of the slices 27r(n — l)i? < t < iirnR, 
crossing the boundary of the slice implies going back from that end of the string to the other and hence 
the situation repeats 
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described by (2.44). So in the compactified time we have the ambiguity that which one 
should describe the fcrmionic wave function in the region < t < 27ri?. To resolve recall 
the macroscopic loop operator W ( 2.15) and the constraint (2.25) are symmetric under 
z+, z+ Z-, Z- . So in the Schrodinger equation (2. 12) we have the symmetry, 

t 2t:R -t; W ^ -W, (3.33) 

which takes an wave function z^,z^ -2- representation^"'"' along with reversal of the 

sign of the energy E — )■ —E in the free fermionic wave function as introduced in (2.28). As 
both z+, z+ and Z- describes same wave function in different representation so it must 
be a symmetry inside the double delta well (note this is never a symmetry outside the 
well where the fermion can see only one potential barrier). The condition (2.25) remains 
unaffected by this symmetry and we can project the wave function t = to = Q = 2'kR to 
the physical sector. Now under the transformation (3.33) the free fermion wave function 
(2.28) is just changed by a phase e^*-^'^^ whereas according to (2.28,2.30,2.41) the wave 
function describes a state \E ±ni) goes from e~*-^*='="* — >■ Qi'E{'^'^R-'t)T{'^'^nR-nt) _ Although 
the transformation (3.33) take the wave function from 2;+ to z^ representation but the 
both have the same time dependent part. So when we consider the wave function (2.41) 
we see it does not respect the symmetry. Note unlike the closed string momentum modes 
which respects the symmetry in Euclidean time because of their periodicity on the circle, 
W{z±z±,Ho) generates discrete shift in energy E ± ni at any R, which leads to the 
violation of symmetry. So we conclude the wave function inside the well which is our 
compact time < t < 2t^R must corresponds to that of a free fermion (2.28). We will 
come to the string theoretical interpretation in the next subsection. The wave function at 
t — Q = 2t^R corresponds to (2.41) so the partition function corresponds to the transition 
amplitude 



lim(V^„(6)|V^>(27ri?-e)) = {^^'^^^^ -^^^-e^^*")]^,) 

= Tr^|.[e-^'^^'^^°[e^(*°)] ] 



(3.34) 



where Tr^B implies the summation over all free fermion eigenfunctions and '0> corresponds 
to (2.41). As /9 — 7> 00 at double scaling limit, so inside the partition function we can replace 
it by '0> e~^^*°^ipo Note when we consider Schrodinger equation, the contribution 



-"^^this is from the definition of z± (2.6) and the reversal of the sign of W is explained from (3. 32), the 
transformation changes the sign in the r.h.s of (3.32) because the time interval e — > — e under the 
transformation and hence the relation will remain unchanged. Note this transformation is also associated 
with the reversal of the sign of the gauge field A. but we chose axial gauge a = a,where a, a are the zero 
modes of gauge fields A, A, so this is not affected 

-'^^In order to reach from the Irst to 2nd step in (3.34)wc utilize the fact that we can scale the time 
t ^ (3t so that the term with the macroscopic loop operator J dtW{t)d{t — to) will get a factor ^. 
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from tlie term ^^^^ W {to) in the expression of hamiltonian (3.34) at double scaling limit 
will not be negligible due to a transformation of variable which leads to the physics at the 
vicinity of the top of the potential, as we discussed section 2.1 and can be found [1]. Prom 
(2.15) we know that in a single fermionic state the presence of D brane implies implies 
the insertion of the following operator 

^wiu) ^ ^ ^ t+i++t-i- + t-i+ + M- _ 

Here first we will evaluate the partition function for single fermionic d.o.f in order to 
understand the behaviour of the system in the presence of a brane. Next we will derive 
the grand canonical partition function. 

Now following the discussion in section 2 we can represent the wave function (2.40) 
Sit t — 2'kR — e in (3.34) with the operator W expressed either in terms of z+, z+ or 
Z-,Z- representation. We have shown in the Appendix that {z+z+\z+z+\z+z+) and 
(^_2;_|f_i_|z_2;_) diverge. So we must express W as l^(f_i_, Hg) for the basis E) 
basis and vice versa. So applying (A. 7) we can write (3.34) as: 



Tr^o 



= Tr, 



^-2nRI3H', 



(3.36) 



Where 

with E'^ is the eigenvalue of iJ^(note that we omitted /3 factor from the expression of 
E'g following the discussion of section 2, which can be done at double scaling limit by 
redefinition of the variable). Before any further analysis let us make the comment that 
here we have considered the macroscopic loop operator in the NS NS sector. For any 
other sector we can use analysis of section 2, expressing W in terms of Z-.Z-{z+, z+) in 
z+, z+{z_, z_) representation and using (A. 7) to express W{to) complete in terms of 
within the trace. Although we will have a very different expression of H'^ but the analysis 
will remain same. Also if we did not apply this condition (2.25), note we will have the 
term (z_i_)_( {z+z+)+ ) in the partition function from the expression of W. This gives 
rise to an infinite contribution to the partition function limr_^oo re"^^^~'^ (where is the 
phase of the wave function) and so the partition function diverge. This is the signature 
of the presence of an unphysical degree of freedom leads to instability of the system due 
to leakage. Now according to the discussion of section 2, at the double scahng limit the 
partition function (3.36) can be expressed as the sum over E'^the eigenvalue of H'^ as: 

.2nRP[E-^log{l-m. 



E 

E 



(3.38) 



Hence in the double scaling limit where ^ ^ oo and with Euclidean time, we can lift up the term 
to the exponential and the exponent gives an exact expression what we have obtained from the path 
integral (3.23) . Also following the discussion of section 2 we can directly add W{to) in the expression of 
hamiltonian in Euclidean time to get the expression (3.34) 
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Now note that E'g,the eigenvalue of H'^, has branch cut at = ^ so we need to subtract 
a small cut-off loge in order to have an well defined expression of the energy and after 



subtraction E ^ E' is an one to one mapping, Also note at the singular point, E = -^f^Bi 
g-27rR/3£;' jg trivially zero and so it will not contribute any pole to integrand. Now the string 
theory compactified at time interval 2t:R is described by the grand canonical partition 
function of fermion at finite temperature and chemical potential ii. So the free energy 
F — logZ in the presence of Dbrane is given by 

[•oo 

J-(At)=/ ci£;p(£;)io^(l + e-^('^+^(^))), (3.39) 

J oo 

wherep(£') is given in (3.17) 

e^ME) ^ R^E) = ^}'^^^^^\ . (3.40) 
^ ' T{-iE + 1/2) ^ ' 

Now we can calculate free energy J-'(/x, R) — logZ(/x, R) as. 

/OO 
dEp{E) log[l + e-^('^+^'(^))] . (3.41) 
-oo 

with the density (3. 17), and from (3.39,3.37) the free energy is given by 



PC 

J — 1 



I g-27rR/3(M+E'(£;))) 

ir 



r=n+|>0 



, 7r* 



r=n+^>0 



where we have 



(t>',{E') = UE) (3.43) 

So from the expression of the density of the energy eigenstates (3.17) (ignore the A factor 
) this imphes the number of energy eigenstates between E to E+dE is same as that of 
between E' to E' + dE' . So the partition function on a circle in the presence of D brane 
corresponds to a deformation in static Fermi sea where the deformation is expressed as 
E — —jj. =^ E' — —jj. with all the energy eigenstates are in one to one mapping. Note 
that the partition function is getting contribution only from the deformation of Fermi 
surface instead of excitation modes. This is because we are in compact dimension and its 
only the trace over energy eigenstates contribute, which we will explain more in the next 
subsection. Finally the expression of the free energy J^(/i, R) in (3.42) along with (3.43) 
suggests the following 

-2^i?/3{i/„-5^1og(l-2«^^)} 



2Ho 
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Tr^,[e-'^^^^''''\ (3.44) 



where 



f(H,) =log{l--^) ; H'^^Ho-,^f{H,). (3.45) 



and in the last step we made a transformation from the basis i^fiE) — e^'^'^°^^^e~'^^*{z±z±) 
il^tiE') = eT^^°(-^')e-^-^'*(^±z±)^^-^'"5 with E' ^ E- ^\og{l - ^) and also (t)'^{E') = 
(f)o{E). Note as we discussed in section 2, although the contribution from the macroscopic 
loop operator has a factor ^ however in the double scaled hamiltonian it cannot be ignored 
and the shifted energy will be given by E'. Above expression implies that the Type OA 
MQM on a circle in the presence of a D brane can be viewed as a free theory with the 
free hamiltonian Ho with the wave function replaced by the above one. This point will 
be relevant in section 5. Note the relation (3.21,3.22) can be expressed as 

^^^^■J'i^^) = €i-^^)- (3-46) 

Also note that from nonlinear relation between the effective hamiltonian H'^ and the free 
hamiltonian Hg its evident that in the genus expansion of free energy in the relation (3.22) 
will have both odd and even powers of /x( this is because the new free energy corresponds 
to replacing // — >■ E{E')\e'=ij. which follows from (3.39)). This is the signature of the 
presence of surface with boundary which is the implication from MQM/string theory 
duality. 

3.4 String theoretical interpretation 

Let us very briefly say about the string theory side of the above story. The matrix model 
partition function (3.44) corresponds to the disk amplitude 

^^|^(L„+r„)|j^ (3.47) 

Where \B) stands for boundary state and the above expression resembles (3.34) on the 
matrix model side. The invariance of theory under the symmetry (3.33) is related to the 
symmetry of the boundary state (2.47,2.48) under interchange of the left and right moving 
tachyonic components in W 

U{k)e'''^^+^^-'^^ O U{-ky^^^-^'^-^'>' (3.48) 

as well as other Virasoro primaries obtained from discrete shift of matter and Liouville 
momentum of tachyonic modes. One can see that in absence of D brane this is a symmetry 
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in MQM. The symmetry is implemented from the interchange X — )■ 2ttR — X followed 
by a sign reversal of momentum k —k which is in accordance (3.33) with k = ^ and 
X — )■ it, U(k) is the wave function. This keeps the matter part ( temporal part in matrix 
model) of the operator/state invariant. The reason for this symmetry in matrix model is 
that the tachyonic states obtained from the 2nd quantized free fermionic theory (ground 
state) or the collective field theory are actually in one to one to one correspondence with 
the above operators(3.48), as we explained in section 2.3. Consequently we found that the 
symmetry (3.33) projects the Hilbert space between <t < 27ri? to free fermionic ground 
state. This is because the discrete tachyonic modes arises from Minkowskian theory with 
shift X ^ iX ; k ^ ik =^ ^ikXf^^ikt^ _^ ^ikx^^^ikty However the states created by the 
action of the operator (^_^+)"(z_z+)'" on ground state, although describes a state with 
imaginary energy but remains in Minkowskian time (unperiodic) and so projected out 
by the symmetry within Q < t < 2t:R. However these states arise on expansion of 
W{z-Z-,'Zj^Zj^,Ho) which actually causes the excitation of a fermion over free fermionic 
ground state,creates a coherent state. 

So we have only free fermionic ground state in the region < i < 2t:R which in string 
theory indicates that we have only free closed string modes in the region < X < 2t^R. 
Coherent states are strongly locahzed at the point of insertion of W. In matrix model 
partition function, the absence of the excitation modes has an explanation in the fact 
that on a circle in the presence of a brane, partition function corresponds to transition 
tp Wip exactly at t = to — = 27rR, where {ip) is the free fermionic ground state. So 
the operators from W which generates excitation naturally will not contribute in the trace. 
This, alongwith the condition (2.25) (which ensures the conservation of fermion number) 
implies that the partition function will corresponds to that of a deformed Fermi surface. 
The poles of the partition function (3.42) corresponds to free closed string tachyons in 
the bulk. In section 5 we will see that same features will be reflected even when we turn 
on the tachyonic background. 



4 Fermionic scattering and semiclassical analysis 

In this section we will study scattering of fcrmions in the presence of D brane and tachyonic 
background at quasi classical limit. The scattering amplitude is given by 

S ^ {/3,t ^ oo\a,t ^ -oo), (4.1) 

where a and (3 denotes the incoming and outgoing state. As the single incoming and 
outgoing state is given by and \z-Z-) respectively [28], so 

S = {z_z_,ovLt\z+z^,m). (4.2) 

Now note that z^z^ and z^z^ representations are related by a unitary operator S, which 
in our case is nothing but the Fourier transformation on the complex plane. Recall, the 
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energy eigenstates in absence of the D brane are given by (2.28). The wave functions in 
(2;+,^+) and {z_,z_) representations are related by 

= j dz+dz+K{z^, z+)K{z^, ^+), (4.3) 

where K{z-, z+) — -|=e*^-^+. Acting on energy eigenstates, we have 



s^^ = 7^(£;)^f , n{E) = (4.4) 



The factor Tl{E) is a pure phase 



n{E)n{E) = n{-E)n{E) = i, (4.5) 

which proves the unitarity of the operator S. Now in absence of the D brane wave 
function (2.28) evolve according to free hamiltonian Hg so from orthonormality of the wave 
functions (4.2) is given by TZ{E)5{E^ — EJ). The operator 5' relates the incoming and 
the outgoing waves and therefore can be interpreted as the fermionic scattering matrix. 
The factor 'R-{E) is identical to the the reflection coefficient. This condition can also be 
expressed as the orthonormality of in and out eigenfunctions 

= 5(E+-E_), (4.6) 

with respect to the scalar product. We usually absorb the factor ??.(£') in phase by defining 

to make the wave function biorthogonal, where (f){E) is the phase of the incoming and 
the outgoing wave function(2.28). Now consider the presence of D brane. For a single 
fermionic state, (4.2) is given by 

{z+z+,t^ oo|e^(*)|z_^_,i = oo) = (^+z+,out|(l- ^"^^^ ^ ~ ^"^° )|^-Z_, in). (4.8) 

Now according to (2.32) W will be expressed either in z+z+ or z-Z- mode 

(^+,£;+,out|e'^(*°)|^_,£;_,in) = {z+, E+\{1 + "-)t=tjz^, EAn) 

= {z^,E^\{l + ^)\z.,E_) 

H'B 

2E 

= {z^,E^\{l-—)\z_,E_) 

H'B 

logil-^-) 

^n(E+)e ^5{E+-E_), (4.9) 
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where in the 2nd step we have out|2;_5_|2;_z_, in) = from (A. 6) . Now 

the presence of the D brane will modify the phase of the outgoing state over the incoming, 

which is given by the factor e '"s . So for the change of phase 5<j){E) we can write 



e 2 = e '^s (4-10) 

Now in (4.10) using the relation (4.7) will leave us with the amplitude e^*^^-^. In [43] its 
explained that the complex phase in the wave function is the signature of tunneling and 
we can presume the above factor accounts for the same. Also note that instead of the 
above macroscopic loop operator in NS sector, if we took the macroscopic loop operator 
in some other sector given by W' according to the discussion in Appendix 
we will have the phase shift of the outgoing state + iW'(l — ^). 

Lets consider scattering of a tachyonic state (which are being created from the action 
of (^±i±)" on fermionic ground state) from D brane. This is better understood from the 
collective field theory where scattering to a single tachyonic state with energy E is given 
by ~ U{E) where U(E) is given by (2.48). This supports the fact that the D brane act 
as a coherent source of closed strings. Now we consider the classical limit /9 — >■ oo. At 
this limit the ground state of MQM is obtained by filling all energy levels up to some 
fixed Fermi energy which we choose to be Ep = —/j,. Quasiclassically every energy level 
corresponds to a certain trajectory in the phase space of variables and they 

are separated by a factor ^. The Fermi sea can be viewed as a stack of all classical 
trajectories with E < Ep and the ground state is completely characterized by the curve 
representing the trajectory of the fermion with highest energy Ep. For the Hamiltonian 
Ho all trajectories are hyperboles Ij^z^ + = —E and the profile of the Fermi sea is 
given by 

z+Z- + 'Z-Z^ = —fi. (4-11) 

First consider the theory without D brane. Then the low lying collective excitations are 
represented by deformations of the Fermi surface, 

z+z_ + z_z+ = M{z+z_ + z_z+). (4.12) 

In order to study the scattering with such deformed background we will follow the analysis 
of [16]. The perturbed wave functions are related to the old ones by a phase factor 

i^i{z±z±) = e=F^^±(^±^±^^Vf (^±^±), (4.13) 

whose asymptotics at large z±z± characterizes the incoming/outgoing tachyon state. We 
split the phase into three terms 

V±{z±z±; E) = y±(z±z±) + ]^(t>{E) + v^{z^z^- E), (4.14) 
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where the potentials V± are fixed smooth functions vanishing at 'z±z± = 0, while the 
term v± vanishing at infinity and the constant (p arc to be determined. Now in order 
to understand the time-dependent profile of Fermi sea first consider the situation in the 
absence of the brane as described in [16]. 

Jo Zj^Zj^\J Z—Z— 

(4.15) 

where M is the normalization. At ^ — > oo Fermi profile can be obtained from saddle point 
approximation which is given by 

'z^z- + z^'z- = -£'± + {z±d± + 'z±d±)i{}±i^±z±). (4-16) 

So following [16] it appears the perturbed state will be an eigenstate of the deformed 
hamiltonian H = Ho + Hp where Hp is given by 

Hp = (^±9± + z±9±)(^±(^±z±; H) (4.17) 

Now in the presence of the D brane scattering matrix element will be given by 

{z+z+, E+, ont\e^^^°^\, E_,z_z_, in) 
. S-matrix element is expressed as 

Sperturb = e-'^AA 7^^£z^^£^e'(^+-+^+-)e-v^(--)(^_^_)*^- 

{ Vz+z+^yz_z_ 

[e-*ito [e-^^p-^(*°)]e*/-oc^°]e-^^+(^+^+) (z+z+)^^+ 







^/z^^/zIzZ 

(4.18) 

So in the presence of the D brane, in the classical regime from (4.18) we can write the 
Fermi profile in the presence of D brane as 

z+z- + z+z- = -E± + {z±d± + z±d±)^±{z±z±) + {z±d± + z±d±)W{z±z±, E) (4.19) 

The perturbed hamiltonian for the deformed state is given by 

H'p^Hp+ {z±d± + z±d±)W{z±z±; H) (4.20) 

So essentially the Dbrane act as a source at to- So as we see in the presence of D brane 
Fermi profile develops instabihty. 
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5 Perturbation by momentum modes 



5.1 Collective field theory analysis 

In this subsection we will consider type OA matrix model with the time t compactified on 
a circle of radius R, perturbed by momentum modes in the presence of D brane. We 
know that the presence of D brane change the tachyonic background [18] so that we need 
to go through collective field theory analysis to understand how the MQM wave function 
in a perturbed background with a D brane is related to the one without a brane. 

Here first we briefiy review the scenario without D brane and then study what happens 
when we consider the theory in the presence of D brane. The tachyon modes of the 
closed string theory are the asymptotic states of collective field theory [25] . The discrete 
tachyonic operator Tn ~ /woridsheet e^*"^/^e*^l"'l/-'^~^^'^ corresponds to the following operator 
in matrix model [26,28], 

V±n/R = e-^'iz^z^r"". (5.1) 

These operators creates a discrete tachyonic state of momenta ^ over the matrix model 
ground state and are periodic in Euclidean time. 

Ti 

[Ho,V±n/R] = Tt^V±n/R \ k > 1. (5.2) 

So V±n/R shift the energy E ^ E^i^ cause a time-dependent perturbation to Fermi sea. 
The perturbed state in general can be expressed as 

*| = eT^'^^^i^'i'-^Vfi = >V±^/^f±, (5.3) 
where the phases (f± have Laurent expansion 

Mz±z±; E) = h{E) + t±k{z±z±f - RJ2 lv±k{z±z±)-''/''.. (5.4) 

^ k>l k>l 

t±k parametrize the asymptotic perturbation by momentum modes of NS-NS scalars, cor- 
responding to the operator introduced (5.1), Note the above wave function asymptotically 
behave as 

^i{z±z±) ~ {z±z±f'''~^ e^^^-^^-^) e'^±{z±z±) . u±{z±z±) = J2 \z±z±\^. (5.5) 

k>l 

Prom the above its evident that tachyonic perturbation can be achieved by deforming the 
integration measures (i[^±^±] to [14] 

[cfl±2;±] — >■ [dz±z±\exp {±iU±{z±z±)) . (5.6) 

Extending the discussions of section 3, these wave functions (5.3) diagonalizes the de- 
formed kernel (5.6). While the perturbed wave function evolves in time with Hg, but 
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it can be seen as the stationary state w.r.t an effective hamiltonian H = Ho + Hp{H), 
where the expressions for perturbed hamihonian Hp from semiclassical analysis is ob- 
tained in [16] as we reviewed in section 4. The partition function is given by Tre'"^^^^^ , 
following section 3 which can also be expressed as Predholm determinant. We have the 



Having given a brief review of type OA MQM perturbed by the momentum modes lets 
consider the theory with D brane. Note that type OA MQM without any net DO brane 
background charge, in the double scaling limit can be viewed as a pair of noninteracting 
fermions moving in inverted harmonic oscillator potential in x and y direction respectively 
(2.8). So from the analysis of [25] its apparent that the respective collective field theory 
will be the generalization of the one for c = 1 case to two (target space ) dimension (x,y) 
and the collective field expressed as 0(x, y, t) = (f){z, z, t), give the eigenvalue density in two 
dimension with appropriate normalization. Collective field hamiltonian will be the sum 
of the hamiltonian for (f){x,t) and (j){y,t). The fiuctuation of the collective field over the 
static value = 0^ + 9^77(t, t) ,where tt^o = Po = Va«f - x'^ - y'^ with pi = pl^ +ply ; p^^ = 
^ ; Poy — ^ and d-,-ri{T, t) — V where '4> corresponds to the 2nd quantized fermionic 
field. Now the presence of D brane implies, inclusion of a macroscopic loop operator 
W (Laplace transformed, locahzed in time) to the collective field theory action which 
essentially creates a coherent state over MQM ground state. In linearized approximation 
we have W{z, to) ^ J dt J dre~''^^dT-r]{T, t)S{t—to). So the presence of the D brane implies 
a field independent source term in the collective field equation of motion (which can be 
viewed as the back-reaction due to the D-brane). Consequently we obtain the solution 
for collective field at t > to as 4) = 0frcc + 0±perturbed where ^perturbed = / drjito, r)G{to, r) 
with j{to) is the current associated with the macroscopic loop operator and G is the 
Green's function. As ^perturbed is independent of so we see that the effect of macroscopic 
loop operator in the action is to change the momentum associated with rj by the external 
current or in other words the classical solution for rj will get a field independent (but 
profile {x{to),y{to)) dependent) shift. The stationary field 0o will also be shifted due 
to the change of the potential. Essentially from string theory side we can just identify 
the current j (transformed from r to space, and Fourier transformed to momentum 
space) with the overlap amplitude {V{p)\B) and the interaction term introduced in the 
collective field action / dtdr j drTjSit — to) as / (j)c\{p){V\B) where 0ci can be viewed as 
closed string field. On this identification we see that quadratic action for t] [25,28] along 
with the source term resembles the closed string field action in the presence of D brane 
S — 'S'ciosed(0ci) + <Pc\{^-,4>)(y{^-,4>)\B) where 0ci is the closed string field. Change in 
momentum of 77 due to interaction with the localized source has an explanation in the 
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fact that closed string momentum is not conserved in the direction of spacehke boundary 
condition of D brane. Now the meaning of the constraint (2.25) is that we have to 
ensure the fact that while the collective field is in interaction with the localized source 
Lint = / drdt j {drri){t.T) 6{t — to), the hamiltonian for 77 will always remain conserved 
which implies time translation invariance of complete action with no leakage from the 
bulk by making StLi^t\t=to = . So the quantized action for rj always gives the propagator 
have a pole corresponding to 2D massless scalar [25] i.e resembles the tachyon. So from 
string theoretical point of view we see that the constraint act as a no leakage condition 
ensures bulk conformal invariance as we mentioned in section 2. 

. Now lets find out the exact form of MQM wave function in the presence of D brane 
in the background perturbed by momentum modes from Collective field theory(which in 
the absence of Dbrane is given by (5.3) ). Lets proceed in the following way. 

The collective field equation of motion implies the classical solution for left and right 
moving field ax± , OLy± (given by (f){x) + a'dx^{x) , (f){y) + a'dyU.{y) where H is the 
collective field momentum), correspond to fermionic momentum density at the edge of 
the Fermi sea [26,28]. Prom the expression of the hamiltonian (2.8), Fermi surface is 
described by 

\{pl+pl)-^,{x^ + f)^-^x. (5.8) 

The collective field equation corresponds to two separate equation for two fermions de- 
scribed by x,y as no interaction exists among them. So the momentum density p±y.-iP±y 
at the edge of the Fermi surface evolves with time t as [26] 

dtP±x + P±xdxP±x - a; = ; dtp±y + p±ydyp±y - y = 0. (5.9) 

If we consider the fiuctuation of collective field around its classical solution PoiPox,Poy)- 
a± ^ Po + e± defining e±. — ^C±) C± is shown to correspond the right and left moving 
tachyonic fiuctuations [28]. Now at classical limit the macroscopic loop operator con- 
tributes a source term to the above 



dtP±x+P±xdxP±x-x 



dW{x,y) 

This essentially gives a discontinuity 



6{t - to) ; dtp±y+p±ydyp±y-y 



dW{x, y) 



dy 

(5.10) 



6{t 



_ dW{x,y) ^ 

P±Xi\to + € — P±Xi\to-€ — \to ) K"^-^^) 

where in (5.11) Xi stands for x,y for i = 1,2. So the effect of D brane is to change the 
Fermi profile above t > to to Po ^ p'g — Po + ^^^^(to). As above to, P± evolves according 
to free hamiltonian Ho so we can express the fluctuation of the collective fleld for t > to as 
e'_j_ = to see ^'^ corresponds to tachyonic fluctuation mode. However the redeflnition 
above to implies a nonlinear shift of collective fleld momenta p'^ = p'o{Po,x,y) and the 
fluctuation = C'{ijPo,x,y). e'^.^ ,e'j.y combined to give complex tachyonic fleld. 
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So from the viewpoint of collective field scenario lets find out the picture in MQM. 
The shift of fiuctuation mode above ^ — > ^' implies a shift in the perturbing phase (5.4) . 

where 

Vw±iz±Z±] E) = 1-4>{E) + t±k{tm±, Vn±, ^, E) ftk{E , Z±Z±) [z±Z±f^^ 

^ k>l ^ 

1 T 

~ -r^±k{tm±,v„±, —, E)fyk{E,z±z±) [z±z±]~'''^ , (5.13) 
k>i ^ 

fyk and /jfcis the extra factor arises due to the nonlinear shift in ^ from z±z± factor 
which arises due to the action of W. These factors give nonperiodic shift to momentum 
by integer numbers i.e | ^ ^ + n with t±k = t±k{tm±,Vn±, E) , v±k = v±k{tm±,Vn±, E). 
Note that once we try to interpret the consequence of (5.11) in MQM we need to replace 
W — )■ Wproj (2.41). The fact that in absence of W we get back our original wave function 
so (5.3) implies that the shifted dressing operator W is of the form W' = F(W)W with 
F — )■ 1 in absence of W. Exploiting the fact that (2.39) is just the quantum version of 
(5.11) and the indication of the semiclassical analysis (4.19) implies W" can be expressed 
as : 

= e^^^'"(^±^±'^)V'f = (1 T W^proj)e^^'^(^±^±'^Vo'L = (1 T W^proj)>VVf±, (5.14) 

This is because at double scahng hmit ^ — > oo the wave function is effectively given 
by e~^yVipg±{as we have seen in (3.34) ). This is the solution from (2.40) if the initial 
free fermionic wave function (2.28) is replaced by the dressed one (5.3). Finally as before 
we will express Wproj in terms of z^z^ {z+z+) in z^z+ {z^z^) basis as this will express 
(5.13) in the following form 

^wp±{z±z±; E) = ^(l){E) + RY,t±k{tm±,Vn±, ^,E)ftk{E,ar{z±z±y'') [z±z±]''^^ 
2 fe>i R 

~ ^Yl TV±k{tm±, Vn±, ^, E)fyk{E, br{z±Z±y) [z±Z±\~''^^ , ' (5.15) 
fe>l 

where r is a positive integer. In next subsection we will see that this expression of </?' lead 
to convergence of the partition function. The tachyonic perturbation can be introduced in 
the path integral by deforming the kernel as in (5.6) and consequently the string partition 
function can be expressed as the Fredholm determinant as in (3.27). We can evaluate the 
Fredholm determinant with a set of diagonalizing wave function which is given by (5.12). 
In the next part of this section we will evaluate the gran canonical partition function in 
the hamiltonian formalism. We will show that the tachyonic deformation in the presence 
of the D brane is generated by a system of commuting fiows Hn associated with the 
coupling constants The associated integrable structure of the partition function is 
that of a constrained Toda Lattice hierarchy. Now in order to see the Toda structure of 
the partition function we need to review Lax formalism. 
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5.2 Lax Formalism 

Here we will briefly review the Lax formalism in the context of Type OA matrix model. 
Consider the operator {z±z±) which can be represented as shift operators o)^^, where cu 
acts on energy eigenstates as cj'^^tp^ = ip^'^\ We have cu = e"""^^ shifts the variable E by 
i. The operators a) and E satisfy the Heisenberg-Weyl commutation relation 

[uj, -E] = ico, [uj-\ -E] = -iu;-^. (5.16) 

Now let us consider the representation of these commutation relations in the perturbed 
theory. The dressing operators >V± (5.5)are now exponents of series in u with ^'-dependent 
coefficient 

The operators 

M+ = -yv+Eyv+^ M_ = -yv-E,yvz\ (5.18) 

known as Lax and Orlov-Schulman operators satisfy the same commutation relations as 
the operators lj and E 

[L+,M+]^iL+ , [L_,M_] = -iL_. (5.19) 

The Lax operators L± represent the canonical coordinates z±z± in the basis of perturbed 
wave functions 

{E\e^'^W±L±\z±z±) = {E\e^'^W±%±z±\z±z±), (5.20) 

while the Orlov-Shulman operators M± represent hamiltonian Hq — —^{z+z^ + z^z^). 
Therefore the L and M operators are related also by 

M+^M_, [L+, L_] = 2iM±, {L+, L_} = 2M| - ^. (5.21) 

The last identity is not satisfied automatically in the Toda system and represent an addi- 
tional constraint analogous to the string equations. The operators M± can be expanded 
as infinite series of the L-operators. Indeed, as they act to the dressed wave functions as 

{E\ e^"^°W± M±\z±z±) = ±i{z±d,^ + z±d,^ + lM{z±z±) 

E ktik{z±z^f'' + + E v±kZ±z^-'A ^|((^±^±)). (5.22) 



^k/R 

,fe>l fe>l 



we can write 

-k/R j^-k/R 



M± = E kt^kL'J'' + /X + E v±kLr . (5-23) 

k>l k>l 
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In order to exploit the Lax equations and the string equations we need the exphcit form 
of the two operators. It follows from that L± can be represented as series of the form 



L_ = e'^'^ i^-^ + a_feC^-i+'^/«j e-"^/2. (5.24) 
Recall that the dressing operators >V± in terms of E and a) are of the for 

= e^"^''' [l + Y,w^k(^^^''^ e^'^^^>i (5.25) 
Studying the evolution laws of the Orlov-Shulman operators, one can find that [20] 

It means that there exists a generating function Ts[t] of all coefficients v±k 

(5.27) 

It is called r-function of Toda hierarchy. It also allows to reproduce the zero mode (p and, 
consequently, the first coefficient in the expansion of the Lax operators 

e^^(^) = r\s - i) = ^^iij^. (5.28) 

We are going to show that the partition function coincides with r-function (5.27). Finally 
note as the partition function is described in terms of the Fermi level /x. So in the 
description of Lax formalism we will replace E by /x. Now let us discuss about the 
integrable flow. Let us identify the integrable flows associated with the coupling constants 

at„L± = [i^„,L±], (5.29) 
where from (5.18), the operators are related to the dressing operators as 

= {dtM^+' = {dtMWz\ (5.30) 

it is clear that Hn — W+uj"'/^W^^+ negative powers of u^^^, which implies expression of 
Hn can be given by [14] 

H^n^{L'^\ + l{L'^\, n>0, (5.31) 
dtrnHn - dt^Hm - [Hm, H^] = 0. (5.32) 
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Equations (5.32,5.30,5.31) imply that the perturbed theory possesses the Toda lattice 
integrable structure. The Toda structure implies an infinite hierarchy of PDE's for the 
coefficients v„ of the dressing operators, the first of which is the Toda equation for the 
phase = (/){£ = — /i) 

^ ^ QiHt^yiHt^-i/R) _ Qi<t>{lJ-+'i/R)-'i<t>{lJ) _ (5.33) 

5.3 String theory on a circle with the D brane in the presence 
of tachyonic background 

In this section we are going to evaluate the free energy of type OA MQM in the presence 
of D brane and with tachyonic background, in the grand canonical ensemble and try 
to understand the relevant string theory. Recall in section 3 we have seen that in the 
absence of the momentum modes, within the time circle < t < 2t:R the solution 
of the Schrodingcr equation corresponds to the free fermionic wave function. This in 
the string theory side giving a picture that we have free closed string states along the 
circle and the coherent states are strongly localized ai t = to = iX° . So with the 
same view in the presence of tachyonic background, within the circle < t < 27ri? 
the wave function must be given by (5.3). The perturbed wave function, while time 
dependent w.r.t the free hamiltonian Ho it is stationary w.r.t an effective hamiltonian H, 
similar as discussed in section 5.1. So lets consider the perturbed MQM with the effective 
hamiltonian H = Ho + Hp{H) in the presence of D brane. First consider the partition 
function (3.7) where now we replace the integration kernels with the deformed measures 
(5.6). So as a generalization of (3.7), in the perturbed background, the Matrix model 
partition function in the presence of D brane (3.23, 3.24,3.25 ) will be with the deformed 
kernel as 

OO jy 

n [dz+,] [dz.,] [dz+,] [dz.,] [citfc]e^*"±(^±^±) « det,, (6^*.-^+.-.) 

— OO 

detjk (e-^«^+^^-^) detjk (/'^'^'+^'-'^) detjk (e-^«^+^-^-*) 

expElog(l + + ^-i^-i + + (534) 

The partition function will be given by Fredholm determinant Det(l + e~'^'^WK)) (3.27) 
where in order to evaluate the determinant we need to choose the basis which diagonalizes 
K (i.e (3.13) with deformed measure as given in (5.34) ) and we evaluate the expectation 
value of W in the same. In order to evaluate free energy in the presence of D brane we 
will proceed in the following way. First consider the scenario without D brane. Recall 
the expression of free energy which is expressed in terms of the phase of wave function 
[14, 15] (which is of the same form of (3.42), expressed in the absence of brane). In a 
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perturbed background the phase (f){E) will be replaced by that of the perturbed wave 
function (5.3) in the expression of free energy [14]. So for the effective hamiltonian H ( 
H — Ho+Hp(H) where Hpin the semiclassical hmit obtained in (4.17) ) of which (5.3) is an 
eigenf unction, the analysis of section 3.3 implies that free energy of the perturbed system 

oo 

in grand canonical ensemble is given by J-" = log-E with Z = ^ e~^^^^^^ {Tre"'^^^^^ } n = 

N=0 

Det (1 + e"^'^^'^'-'^"'"^^). This is supported from the view of [16] where in the semiclassical 
regime the explicit expression of Z is obtained in this form. The Fredholm determinant 
(3.13) in a perturbed background is given by Z [14]. In the presence of D brane we have 
the Fredholm determinant (3.27) which in perturbed background is expressed in (5.34). 
So as in section 3.3 free energy must be obtained from the thermal partition function in 
the presence of D brane i.e by insertion of the operator e'^*^*"^ in the partition function and 
evaluating the expectation value. So here in hamiltonian formalism we will evaluate the 
grand canonical partition function Det(l + e"^*^^'*"^-') in the basis (5.3) with the insertion 
of the operator and it must be same as the Fredholm determinant (5.34). Here we are 
going to show that if we consider the prtojected theory as described in section 2, the 
above mentioned grand canonical partition function have the integrable structure of tau 
function of Toda hierarchy. Now the partition function with the momentum modes in the 
presence of the D brane is given by the transition amplitude from the initial state Wipo to 
the final stateVV''?/'o>5 where ipoy is given in (2.41) and they represent the fcrmionic wave 
function before and after being scattered from the D brane and the corresponding dressing 
operator is W .Now note that in a compact dimension just before being scattered, the 
wave function at t — 2t:R — e must be given by the one &X, t — e with a time evolution 
2t^R. This leads to the identity 

WMo) = W'(l - W^(t±i±, Ho))MQ = (l - W{%±z±, Ho)) WMto), (5.35) 

where to = = 2tiR. The above relation can also be viewed from (2.40) in the presence 
of tachyonic background, if we replace the initial fermionic wave function (2.28) by the 
dressed one (5.3). Hence the partition function on the circle corresponds to the transition 
amplitude 

Z = lim (1^0(6) |>VVo>(27ri? - e)) 

e— >0 

= lim (W^o(e)l(l - iy(t±5±, i/o))| Wo(27ri? - e)) 

= Trw^„{e-2^^'^^e^(*=°)}. (5.36) 
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Where the partition function is evaluated in Euclidean time and Tr wi/)o denotes the trace 
taken w.r.t (5.3)^^. Grand canonical Partition function will be given by the following 
expression where we will have the contribution from singlet states only 

Jill + e-2^-^('^+^)(V^,^|e^(*=°)|V^^^)} = + e-'^'^^^^+^^i^^le'^li^^)}, (5.37) 

E E 

Note, we could write the above expression for grand canonical partition function only 
because W can be expressed as the direct product of the operators for the single fermionic 
states. Now following (2.15, 2.41) the partition function (5.36) can be expressed as 

n [1 + e-'^^^^^+^\i^f\e'^^^^\i^f)] 

E 

= + e-'^-^(^^^\i^^\{l + '-^^^^^^M)]. (5.38) 

E I^B 

Now we have shown in the appendix that {z^z^\'z^z^\'z^z^) and {'z_z_\'z_z_\'z_z_) di- 
verge. So we must express W as W{z_z^, Ho) for the basis \z^z^, E) basis and vice versa. 
Now note according to the commutation relation (2.30) and from the form of the wave 
function (5.3) 

z_z_ = e 2 -7^^-^ — e 2 

. So according to the analysis Appendix, (V^^l \ i^p) — except when R is an integer^^. 
However for R an integer it contributes a constant term independent of 0, E, in the 
partition function and can be ignored by subtracting out an overall constant from the 
hamiltonian which amounts to multiplying the partition function by an overall factor. 
For the macroscopic loop operator in any other sector, we can proceed in the same way . 
So we can write the partition function as 

n[i + + ^±^^^)e-«''+^))iv^f )] 

E I^B 

= W + (^p'^|e'°'^'"'S^e-2«'^+^))|^f )] 

E 

1 Ho 



E 



E 



Tr^|;[e-2«'^+^°+^-)] (5.39) 



-'^^In order to reach from the 2nd to 3rd step in (5.36)we utihze the fact that we can scale the time 
t /3t so that the term with the macroscopic loop operator J dtW{t)6{t — to) will get a factor so 
that in the double scaling limit where /3 — > oo and with Euclidean time, we can lift up the term to the 
exponential and the exponent gives an exact expression what we have obtained from the path integral 
(2.36) 

^^This is because we can write the integral as {ipo\{'z+z^)^~^\il>o) which following the analysis of 
Appendix- A contributes only at pole. 
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where H'^ is as discussed in (3.45), given by H'^ = i/o~ 2^7^1og(l — ^^); Hp = Hp{z±z±, H) 
is the effective perturbation in the presence of momentum modes 

First note that the eigenstate oi H = Ho + Hp does not diagonalize the complete 
effective hamiltonian Hes — H — 2:??^log(l — 2^). This is the indication from collective 
field theory as we discussed in section 5.1 that the effect of putting a macroscopic loop 
operator in MQM action is to deform the static Fermi sea as well as the tachyonic back- 
ground which shows up as a nonlinear shift of the perturbing phase (p — > (fwp (5.15,5.14). 
Clearly the eigenfunction which diagonahzes the complete effective hamiltonian H'^ + Hp 
will be given by a shift W — > W" = e'^ where i0' will be of similar form of (5.15) and so we 
will denote it by Now we can express W' as the product W' — U{z±z±, Ho)W where 
U is the factor responsible for the presence of macroscopic loop operator so that in absence 
of the operator we should have W = W. Now in order to diagonalize lets recall the tricks 
we used in (3.44). If WVo diagonahze the deformed hamiltonian in (5.39) then we can re- 
place the partition function with deformed hamiltonian H — Hj + Hp{H ,z±z±) evaluated 
in the basis >Wo(-E') with the one in the shifted basis >VVo(-E') Wi/joiE') evaluated 
w.r.t the hamiltonian H — Ho + Hp{H,z±z±). By this shift we can identify U{z±z±, Hg) 
with an operator which has an effect to shift E' ^ E in the eigenfunction ■0^(£"). So the 
operator U will be of the form U{J2n 0'n{z±z±y^ ..)Ui, where Ui is an operator making uni- 
tary transformation to W. Clearly W and W are not related by unitary transformation 
and so Wipo and Wtpo defines the basis of completely different Hilbert space. So in the 
presence of the Brane we need to evaluate the partition function in the basis WVo which 
we can view as deformed tachyonic background caused by the presence of the brane. One 
can verify the expectation value of macroscopic loop operator W{z±z±, Ho) in this basis 
effectively gets contribution from its hamiltonian part i.e W{Ho) as in the case for unde- 
formed basis (5.3) and following similar steps we will get the partition function (5.39) with 
a shifted basis Wipo VVVo- The perturbing phase (f)' of the shifted basis in principle 
can have a complex part for which we need to choose appropriate normalization. However 
in order to evaluate the partition function we will follow (3.44). This partition function 
is exactly given by the one with a shift ^p'^f = W^poiE) =^ tpp' = WtpoiE') and the 
perturbing phase (j)wp{E) — )> (j){E{E')) = (f)'{E') evahiatcd w.r.t the effective hamiltonian 
H = Ho + Hp{H,z±z±) but without insertion of the macroscopic loop operator W, where 
(()wp{E) is as introduced in (5.15) and ijj'^ is the basis which diagonahzes the deformed 
hamiltonian H'^ + Hp^aswe discussed. So following (3.44) we can evaluate the partition 
function (5. 39) in the shifted basis 

^^In order to reach from 2nd to 3rd step we used the same tricks of section 3 which imphes that around 
a delta function in time, we can make the time interval infinitesimally small so that we can ignore the 
commutator terms {[H'g,Ho] +.... higher commutators) what can arise on exponential as a consequence 
of Baker Hausdorff formula 
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where the shifted dressing operator is given by 

and 0(-E') is the pase for perturbed wave function (5.4). So following (3.42) we have the 
free energy given by 

:F{i,,R) = <p{E{E'^ — -f,)) 

= E '^'(^-Z^)' (5-41) 

where we have 

cP'iE') = 4>{E) (5.42) 

So we see the partition function in the presence of D brane in a background perturbed 
by momentum modes with compactified time is obtained from the one without D brane 
by the shift 

E^E' ■ W ^Ws 
which defines a deformed Fermi surface.. 

5.4 Lax formalism for Type OA MQM in the presence of D 
brane 

Our lesson from the previous discussion is that Toda structure for Type OA MQM per- 
turbed by tachyonic modes ,in the presence of D brane can be obtained when we replace 

1 9 F 

V'f ^ = ME - ^ log(l - -^)), (5.43) 

M'+ = Ws+EWs+-\ M_ = Ws-EWs- (5.44) 
Note the operator algebra (5.19) remains same. 

{E\e^"^'Ws±L'±\z±z±) = {E\e^'^'^Ws±^±z±\z±z±), (5.45) 
where 0' = 0(£"). Expression of M' is as described in (5.21) 

{E\ e^"^'W± M'±\z±z±) = ±t{z±d,^ + z±d,^ + 1)^^ {z±z±) 

E k t±k {z±z^f + E' + Y: v'ik {z±z^r'^A *f ((^±^±)). (5.46) 

^,fc>l k>l J 
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As the partition function described in terms of Fermi level /j, 



k>l k>l 



, -k/R 



(5.47) 



The structure of the integrable flow remain same. The Toda flow equation will be given 
by ^ cl>{E' = 

. d d 



(5.48) 



^1 dt-i 

Now in order to see that partition function is a tau function of Toda lattice hierarchy first 

1 



note that 



n>0 



ip(p 1 1 — — — n 



(5.49) 



(5 R 

with Fermi level E' = — /x. Now on the other hand, the zero mode of the perturbing 
phase is actually equal to the zero mode of the dressing operators (5.43). Hence it is 
expressed through the r-function as in (5.28). Since the shift in the discrete parameter n 
is equivalent to an imaginary shift of the chemical potential /x, so (5.49) implies 



(5.50) 



However from (5.28) we have 



2RI3> 



So one concludes that 



'^o(/^ 2Rp) 



(5.51) 



(5.52) 



5.5 Representation in terms of a bosonic field 

Here we will study the classical limit following the analysis of [14] The momentum modes 
can be described as the oscillator modes of a bosonic field ip{z+z+,z_z_) = ip^{z^z^) + 
(p-{z-Z-). The bosonization formula is 

*f ="''"'(z±z±) = Z-ie±^'^±(^±^±) • Z. (5.53) 

(Note here in the presence of FZZT brane n corresponds to the deformed Fermi surface) 
where Z is the partition function and 

Mz±z±) = +i? E tk(z±z±)''/'' + ^d^ + l^ log^i^i - i? E t(^±^±)"'^'' (5.54) 
k>i k>i "'^k 

Then from (5.46) the operators M± are represented by the currents z±z±d±i^ 

Mi-^l{z^z^)\E=-^-i = Z-H^z^d^ip ■ Z. (5.55) 
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5.6 The dispersionless (quasiclassical) limit 

We consider the quasiclassical limit /3 — )■ oo. In this limit the integrable structure de- 
scribed above reduces to the dispersionless Toda hierarchy where the operators /i and Cj 
can be considered as a pair of classical canonical variables with Poisson bracket 

{cu, ii) ~ ui (5.56) 

Similarly, all operators become c-functions of these variables. The Lax operators can be 
identified with the classical phase space coordinates 'z±z±^ which satisfy 

{z+,z.} = {z+,z.}^l (5.57) 

The shape of the Fermi sea is determined by the classical trajectory corresponding to the 
Fermi level E' — —ji. So we have 

-z^z. + z^z. - log(l - ^^^+"- + "+^-^ ) - log 6 = (5.58) 

Where loge is the cut-off cancelling the singular contribution from the point (1 — 
2{z+z-+z+z-) \ _ g ^j^g perturbed theory the classical trajectories are of the form 

Mb 

z±z± = L' ±{u!, n). (5.59) 
where the functions L± are of the form 

L'±(a;, Pi) = e^^''^' u^' (^l + Y. a'±M . (5.60) 

The flows Hn become Hamiltonians for the evolution along the 'times' t„. The uni- 
tary operators VV± becomes a pair of canonical transformations between the variables 
uj.fi and L±,M±. Their generating functions are given by the expectation values S± — 
(p±{z±z±) ■ Z of the chiral components of the bosonic field 

S^^±RY^t^k (^±^±)''"' + /ilog(^±^±)-(?>'±i?Er4 (^±^±)"'^^ (5.61) 
fe>l fe>l 

where V]^ — dF/dtk. The differential of the function S± is 

dS± = M±d\og{z±z±) + logo; d// i? ^ Hndtn- (5.62) 

If we consider the coordinate a; as a function of either ^+2;+ or then 

^ ^ ^d^s+[z+z+) ^ ^d^s.{z.z.)_^ ^5_g3^ 

The classical string equation 



(5.64) 
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can be written as 



Zj^Z^ + Z-^Z^ — 




(5.65) 



k>\ fe>l 



6 Conclusion 

Here we have studied Type OA matrix model in the presence of spacehke D brane which 
are locahzed in matter direction. In matrix model this is expressed by insertion of an 
operator e^*^*"^ into the path integral. When we studied the respective MQM we found by 
application of Ward identity that the time translation invariance of the path integral in 
the presence of such operator gives the signal of leakage of MQM hamiltonian. However 
in dual string theory this phenomenon has a meaning that closed string hamiltonian is 
undergoing a leakage when the string is getting scattered from Dbrane! in order to obtain 
right string theory picture we impose a constraint (2.25) on matrix model path integral 
in the presence of D brane. We explained that this condition has an effect to constrain 
the Hilbert space generated by macroscopic loop operator while keeping type OA MQM 
unaffected. We have shown that when we impose the constraint we get the matter one 
point function from collective field theory. We have further shown that exactly at the 
point of insertion of the brane ( which in string theory correspond to the point where 
open string ends are localized) the wave function for the right and left moving component 
of boundary state with any momentum appears to be identical which can be seen in matrix 
model as a consequence of this constraint. We also found right transition amplitude from 
a free fermionic state to coherent state. Next we consider type OA MQM with the time 
t compactified on a circle. We have shown matrix model path integral in the presence of 
Dbrane can be expressed as Predholm determinant. We evaluated the thermal partition 
function in grand canonical ensemble. As the theory is defined on a circle so the partition 
function correspond to that of a deformed Fermi surface. We have further shown that in 
absence of any such constraint, the partition function diverges. Theory on the circle also 
posses a symmetry (3.33), which is parallel to string theory as we discussed in section 
2.4. This symmetry clearly indicates that coherent states are strongly localized at the 
point of insertion of macroscopic loop operator. Finally we considered type OA MQM 
in the background of momentum modes. First in section 4 we made a semiclassical 
analysis, studied fermionic scattering in the presence of D brane. We found the effective 
hamiltonian in the perturbed background from semiclassical analysis. Its known that the 
presence of D brane change the tachyonic background. So from collective field theory 
analysis we found the right expression of MQM wave function i.e exact modification of 
the dressing operator in the presence of Dbrane. Wc derived the grand canonical partition 
function in the perturbed background in the presence of D brane. We have shown the 
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partition function corresponds to tau function of Toda hierarchy. We have also analyzed 
the theory in dispersionless limit. 

Its interesting to study T duality between type OA and type OB MQM in the presence 
of D brane. One can also study the theory in the presence of flux background and see the 
consequence of the constraint. 
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Appendix 

A Appendix 

Here we will show that type OA MQM wave functions in the presence of D brane satisfy 
orthogonality and biorthogonality conditions. 

A.l Orthogonality Condition 

The wave function is expected to show orthonormal property. 

1. For t < to we have the wave function t/jg given in (2.28). One can check the orthonor- 
mality property by considering the contour integral ,and it is given by [15] 

{iPe'\^e)^S{E-E'). (A.l) 
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For t > to, first consider the wave function ■i/j^{z+, z+,t; E) given in (2.42). we have 

J Ud 



[{1 _ iog(i + ^^^-^-'^^°^ )}e-'^(*-*°)e-^^(^)(z+z+)'^-^]. (A.2) 

Now in order to see the orthonormal property first recall the commutation relation (2.30) 
Note that i+ and i_ shifts E by -i and +i respectively.So we can write 

z\z^ = e^e-^-e^ (A.3) 

Similarly 

z-z-^e 2 e"=^^e 2 . (A.4) 
Now to evaluate (A.2) first consider the expression 

(V'+|[t-i-]"^|V'+) = / d^+rfl+{e^^'*et^°(^')(z+z+)-^^'-^} 
oz+dz+ 

For E ^ E' this can just be written as 

(V'+|[t-i-]"*|V'+) = e-^'^°e'"^^^e5'^'' j dz+dz+{e'^''e-^'f"'^^'\z+z+y''''--^} 

= e-t^°e*'"^-e^^°(V^+(^')lV'+(^))w 

= 0. (A.6) 
For E^E' (A.2) takes the form 

= {e^'^°(-^)e-^^°(^-'"^)} y dz+d^+[z+z+]-"^-^e-'"*. (A.7) 

From contour integral which is for m > 1. Also we conclude that (■?/'+ 1 [^+^_|_]™|-?/'_|_) and 
{ip~\['z^z_]"^\^_) diverge. Here before going to show the orthogonality lets consider the 
situation when m is not an integer. This kind of integration we had in the section 5 in the 
expression {ipoliz+z^) \i^o)- We can consider this as the product over two branh cut 
integrals z^ and z^ and the each branch cut integral can be expressed as the sum of two 
standard contour integral with the cut on right and left side of the real axis > a; > cxd 
and —00 > X > respectively and with a pole at a: = 0. One can see that the integral 
turns out to be zero for any nonintcgcr ^. We have nonzero contribution only when R is 
an integer and the contributing term is m=R.This term corresponds to a pure pole and 
give a constant contribution to the integral, Now back to the question of orthogonality. 
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So using (A. 4 ,A.7), we can be written (A. 2) as 

{,l;iE',tmE,t)) = [ dz+dz4{l - log(l + l^"^ )}e^^'(^-*°)e^'^°(^')(z+z+)-'^'-^] 
[{1 _ iog(l + Z(2l^)}e-^(*-*<')e-^^°(^)(^+^+)^^-^] 
= / dz+dz4{l - log(l + :^i^^)}e'^'(*-*<')et^°(^')(^+^+)-^^'-^] 



Now write 0^ = ^^^je + «0o/m- Shifting ^^/^ ^- - |log[l - log(l + ^^^-^^ ' )] we get the 



lilting, yjolm ~r H'oLm ~ 2^'-'&L-^ ~ ^"-"BV"^ "i 

orthogonal property (A.l) with 

i)^{E,t) = e^^(*-*°)e-^'^°+(^)(z+z+)'''-^ (A.9) 

where 



</.„+(E) = ^E) - zlog[l - log{\ + ^^%^)]. (A.IO) 
Similarly for the wave function in z- representation we find 

<^o-{E) = <^,{E) + %log{\ - log{\ + ^^^%^)] (A.ll) 

So we see the consequence of the insertion of macroscopic loop operator is that the phase 
of the wave function develops an imaginary part which is associated with tunneling. 

A. 2 Biorthogonality Relation 

The wave function is expected to satisfy the following biorthogonality condition which 
has a consequence in the evaluation of scattering amplitude [14-16]. 

y dz^dz-dz^dzJ\)l{:z^z^, t)e''^'+'-+^-'+^iP^' {z_z_,t) = 6{E - E') (A.12) 

Now for t < to we have ip = ipo and for which biorthogonality relation is already derived 
in [14], [15] giving e^'^o(^) = p[^g+|/^j [15]. For t>to, .Biorthogonality relation takes the 
form 

j dz+dz-dz+dz- ^+(^+2+, - W{z+,%+, Ho,t)}e'^'+'-+'+'-^){l - W{z+,%+, Ho,t))ip^' 
= S{E-E'), (A.13) 

In order to show this note that 

J dz+dz-dz+dz- i4{z+,z+,t){l - W{z-,%-, Ho,t)}e'^'+'-^^+^-^){l - W{z-,%-, Ho,t))i/j^' 
= I rfz+rf^+^(;.+ ,^+,t){l-iy(i_,f_,if«,t)}(l-iy(^,^,if„t)) 
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j dz-dz-e'^'+'-+'+'-^)iP^'{z-,Z-,t) 
^ R{E) j dz+dz+i4{z-,z^,t){l-W{z-,Z-Ho,t)}{l-W{z-,Z-,Ho,t)) 

ilj^'(z+,z+,t) 

= R{E)e'^°+5{E - E'). (A. 14) 

Where in order to come from 2nd to 3rd step we used the fact that in 2;+ representation 
we have = ^ The integral in the 4th step, we have evaluated in (A. 8) leads 

to the last step. In order to get the biorthogonality relation(A.13) we must need to set 
gi</.o+(B) _ Compared to the t <to case note the shift of (t)o{E) due to the insertion 

of macroscopic loop operator. 

/OO 
dz-dz-e'^'+'-+'+^-^)ilj^'{z-,Z-,t) (A. 15) 

-00 

R(E) getting absorbed to decide 0o and shifting 0o according to (A. 11) we get the above 
result. 
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